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Abstract 

We explore the parameter space of the minimal supersymmetric extension of the Standard 
Model (MSSM), allowing the soft supersymmetry-breaking masses of the Higgs multiplets, 
mi^2, to be non- universal (NUHM). Compared with the constrained MSSM (CMSSM) in 
which mi 2 are required to be equal to the soft supersymmetry-breaking masses mg of the 
squark and slepton masses, the Higgs mixing parameter and the pseudoscalar Higgs mass 
m^, which are calculated in the CMSSM, are free in the NUHM model. We incorporate ac- 
celerator and dark matter constraints in determining allowed regions of the (/i, m^), (/i, M2) 
and (mi/2,'^0) planes for selected choices of the other NUHM parameters. In the examples 
studied, we find that the LSP mass cannot be reduced far below its limit in the CMSSM, 
whereas may be as small as allowed by LEP for large tan (3. We present in Appendices 
details of the calculations of neutralino-slepton, chargino-slepton and neutralino-sneutrino 
coannihilation needed in our exploration of the NUHM. 



CERN-TH/2002-238 
October 2002 



1 Introduction 



The hierarchy of mass scales in physics is preserved in a natural way if supersymmetric 
particles weigh less than about a TeV. Many supersymmetric models conserve the quantity 
R = (^—iy>B+L+2S^ where B is the baryon number, L the lepton number and S the spin. 
If this R parity is conserved, the lightest supersymmetric particle (LSP) is expected to be 
absolutely stable. The most plausible candidate for the LSP is the lightest neutralino X) 
which is a good candidate [1] for the cold dark matter (CDM) that is thought to dominate 
over baryonic and hot dark matter. 

In this paper, we refine and extend the many previous calculations of the relic LSP 
density in the framework of the minimal supersymmetric extension of the Standard Model 
(MSSM). In particular, we expand the recent analysis of the MSSM parameter space in [2], 
where we allowed non-universal input soft supersymmetry-breaking scalar masses for the 
Higgs multiplets. Here we explore in more detail the constraints imposed by accelerator 
experiments - including searches at LEP, b ^ s'j and — 2 - and the cosmological bound 
on the LSP relic density. 

Before discussing our calculations in more detail, we first review the range of the relic LSP 
density that we prefer in our calculations. An important new constraint on this is provided 
by data on the cosmic microwave background (CMB), which have recently been used to 
obtained the following preferred 95% confidence range: flcDuh^ — 0.12 ± 0.04 [3]. Values 
much smaller than ^cdm^^ = 0.10 seem to be disfavoured by earlier analyses of structure 
formation in the CDM framework, so we restrict our attention to flcBuh'^ > 0.1. However, 
one should note that the LSP may not constitute all the CDM, in which case ^lsp could 
be reduced below this value. On the upper side, we prefer to remain very conservative, in 
particular because the upper limit on ili^sp sets the upper limit for the sparticle mass scale. 
In this paper, we use Ocdm^^ < 0.3, while being aware that the lower part of this range 
currently appears the most plausible. 

The parameter space of the MSSM with non-universal soft supersymmetry-breaking 
masses for the two Higgs multiplets has two additional dimensions, beyond those in the 
constrained MSSM (CMSSM), in which all the soft supersymmetry-breaking scalar masses 
mo are assumed to be universal. In the CMSSM, the underlying parameters may be taken 
as mo, the soft supersymmetry-breaking gaugino mass mi/2 that is also assumed to be uni- 
versal, the trilinear supersymmetry-breaking parameters Aq that we set to zero at the GUT 
scale in this paper, the ratio tan (3 of Higgs vacuum expectation values, the Higgs superpo- 
tential coupling /j, and the pseudoscalar Higgs boson mass tua- Two relations between these 
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parameters follow from the electroweak symmetry-breaking vacuum conditions, which are 
normally used in the CMSSM to fix the values of /i (up to a sign ambiguity) and ttia in 
terms of the other parameters (mo, mi/2, ^o, tan (3). 

In the more general MSSM with non- universal Higgs masses (NUHM), the parame- 
ters II and niA become independent again [4, 5, 6]. Thus one may use the parameters 
(mo, mi/2, A*, TTiA, Aq, tan (3) to parametrize this more general NUHM. The underlying theory 
is likely to specify the non-universalities of the Higgs masses: m^ = sign(m?)|mi/mo| : 

1 — 1,2, so it is important to know how the different values of the rrij map into the 
(mo, mi/2, A*, ^0, tan/3) parameter space, a point we discuss in Section 2. Furthermore, 
this non- universality leads to new coannihilation processes becoming important, which are 
discussed in Section 3. 

We review and update in Section 4 the experimental and phenomenological constraints 
on the MSSM parameter space that we use, applying them to the CMSSM. Then, in Section 
5, we explore the NUHM parameter space. Previously, we gave priority to a first scan of 
the extra dimensions of the parameter space, and postponed a complete discussion of the 
NUHM at large tan/3. Here we also show how our results in the (/x, m^), (//, M2) and 
(mi/2,mo) planes for tan/3 = 10 change at larger tan/3, concentrating on the behaviour of 
the rehc LSP density, but also incorporating constraints on the NUHM from accelerators. In 
our discussions of these planes, we emphasize the novel features not present in the CMSSM, 
such as the forms of the regions in which the LSP is charged, e.g., the hghter r, regions where 
the LSP is a sneutrino u (thus necessitating the inclusion of additional coannihilation 
processes), and the potential importance of Higgsino coannihilation processes. These are 
not usually relevant in the CMSSM, where the relic LSP is usually mainly a B. Section 6 
summarizes some conclusions from our analysis, including comments on the range of LSP 
and pseudoscalar Higgs masses allowed in the NUHM. 

The Appendices provide the information needed to reproduce our calculations of coanni- 
hilation processes relevant to this NUHM analysis. In particular. Appendix A hsts the MSSM 
couplings we use. Appendix B extends previous results on neutralino-slepton coannihilation 
to include left-right (L-R) mixing. Appendix C discusses chargino-slepton coannihilation 
processes, and Appendix D concerns neutralino-sneutrino coannihilation. 

2 Vacuum Conditions for Non-Universal Higgs Masses 

We assume that the soft supersymmetry-breaking parameters are specified at some large 
input scale Mx, that may be identified with the supergravity or grand unification scale. 



2 



The low levels of flavour-changing neutral interactions provide good reasons to think that 
sparticles with the same Standard Model quantum numbers have universal soft scalar masses, 
e.g., for the cl, Jj-l and tl- Specific grand unification models may equate the soft scalar 
masses of matter sparticles with different Standard Model quantum numbers, e.g., (d, sJ))l 
and (e, /J, f)^ in SU(5), and all the Standard Model matter sparticles in SO(IO). However, 
there are no particularly good reasons to expect that the soft supersymmetry-breaking scalar 
masses of the Higgs multiplets should be equal to those of the matter sparticles. This is, 
however, the assumption made in the CMSSM, which we relax in the NUHM studied here ^. 

One of the attractive features of the CMSSM is that it provides a mechanism for generat- 
ing electroweak symmetry breaking via the running of the effective Higgs masses-squared m\ 
and m\ from Mx down to low energies. We use this mechanism also in the NUHM, which 
enables us to relate m\{Mx) and m|(Mx) to the Higgs supermultiplet mixing parameter 
and the pseudoscalar Higgs mass ruA- Therefore, we can and do choose as our indepen- 
dent parameters iJi{mz) = A* and mA{Q) = friA, where Q = {mj^mj^y^^; as well as the 
CMSSM parameters (mo(Mx), mi/2(Mx), Aq, tan/3). In fact, in this paper we set ^40 = 
for definiteness. 

The electroweak symmetry breaking conditions may be written in the form: 



where and A^-^'^^ are loop corrections [8, 9, 10] and mi,2 = mi^2{mz). We incorporate 
the known radiative corrections [8, 11, 12] ci,C2 and relating the values of the NUHM 
parameters at Q to their values at m^: 



m\{Q) = ml{Q) + ml{Q) + 2ij.\Q) + A^(g) 



(1) 



and 



ml - ml tan^ (3 + |m|(l - tan^ (3) + A^ 
tan2/3- 1 + Aif^ 



(2) 



ml{Q) 
ml{Q) 




(3) 



Solving for ml and m|, one has 



mi(l + tan2/3) = m^(Q) tan^ /3 - //^(tan^:^ + 1 - Ajf)) - (ci + C2 + 2c^) tan^ /? 
-^a{Q) tan^ 13 - \mlil - iw? (3) - A^ 



(4) 



For models with non- universality also in the sfermion masses, see [4, 7]. 
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and 



mliQ) - //'(tan^ + 1 + A^f)) - (ci + C2 + 2c^) 
-A4g) + ^m^(l-tan^/3) + AW, 



(5) 



which we use to perform our numerical calculations. 

It can be seen from (4) and (5) that, if ttia is too small or is too large, then ml and/or 
m| can become negative and large. This could lead to m^(Mx) + ix^{Mx) < and/or 
m\{Mx) + n'^iMx) < 0. thus triggering clcctroweak symmetry breaking at the GUT scale. 
The requirement that electroweak symmetry breaking occurs far below the GUT scale forces 
us to impose the conditions m^(Mx) + /i(Mx),m|(Mx) + ii{Mx) > as extra constraints, 
which we call the GUT stability constraint ^. 

Specific models for the origin of supersymmetry breaking should be able to predict the 
amounts by which universality is violated in m^g, which can be read off immediately from 
(4, 5). Alternatively, for a given amount of universality breaking, these equations may easily 
be inverted to yield the corresponding values of /i and ttia- In this paper, we plot quantities 
in terms of fi and m^. 

In the CMSSM, to obtain a consistent low energy model given GUT scale inputs, we 
must run down the full set of renormalization group equations (RGE's) from the GUT scale 
and use the electroweak symmetry breaking constraints which fix /j, and uia- Consistency 
requires the RGE's to be run back up to the GUT scale, where the input parameters are 
reset and the RGE's are run back down. Many models require running this cycle about 3 
times, though in some cases convergence may be much slower, particularly at large tan/3. 
Indeed, there are no solutions for < when tan/5 is large (> 40) because of diverging 
Yukawas. In the NUHM case considered here, we have boundary conditions at both the 
GUT and low-energy scales. Once again, the numerical calculations of the RGE's must be 
iterated until they converge. However, in this case, it is not always possible to arrive at 
a solution, especially for large tan p. In our subsequent calculations, we start by making 
guesses for the values of mi,2(Mx) for use in the first run from Mx down to ruz, and it can 
happen that the iteration pushes the solution away from the convergence point instead of 
towards it. Therefore, the first few iterations must be monitored for any potential blow-ups. 
^For a different point of view, however, see [13]. 
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3 Renormalization and Coannihilations in the NUHM 
Model 



The RGEs for the NUHM have additional terms beyond those appearing in the CMSSM, 
and the resulting sparticle spectrum may exhibit some novel features, as we now discuss. 

The new terms in the RGEs which vanish in the CMSSM involve the following combina- 
tion of soft supersymmetry-breaking parameters [12] : 



S = 4(ml-ml + 2{mj 



m~ — 2r7i~ + m~ + mr- ) 

Ll "i{ dR en' 



(6) 



(rr^ , _ , _ 

QzL LzL tR bR 

Here Ql, Ll are the first two generations left-handed sfermions, and Qsl, L^l are the third- 
generation sfermions. These new terms appear as follows in the RGEs for the NUHM: 
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where the Mi^2,3 are gaugino masses that we assume to be universal at the GUT scale. 

In the CMSSM, with all scalar masses set equal to mo at the GUT scale, S — initially 
and remains zero at any scale [14], since 5" = is a fixed point of the RGEs at the one-loop 
level. However, in the NUHM, with mi 7^ m2, as seen in Eq. (6) 5 7^ and can cause the 
low-energy NUHM spectrum to differ significantly from that in the CMSSM. For example, 
if < the left-handed slepton can be lighter than the right-handed one. Also, mf and ml 
appear in the Yukawa parts of the RGEs for the third generation (Eq. (7)), so NUHM initial 
conditions may cause their spectrum to differ from that in the CMSSM. In the NUHM case, 
depending on the parameters, we may find the LSP to be either (i) the lightest neutralino 
X, (ii) the lighter stau fi, (iii) the right-handed selectron cr and smuon JIr ^, (iv) the left- 
handed selectron and smuon p,L, (v) the electron and muon sneutrinos Ug^^, (vi) the tau 
sneutrino z/,-, or (vii) one of the squarks, especially the stop and the sbottom ^. Note that 
in the cases that we consider here, the z5g^ are generally lighter than the Ur in the regions 
in which they are the LSP (that is, when mf < 0, cf. Eq. (7) and Fig. 2 of ref. [2]), unless 
niA is very large (^ 1000 GeV). 

We assume that R parity is conserved, so that the LSP is stable and is present in the 
Universe today as a relic from the Big Bang. Searches for anomalous heavy isotopes tell 
us that the dark matter should be weakly-interacting and neutral, and therefore eliminate 
all but the neutralino and the sneutrinos as possible LSPs. LEP and direct dark-matter 
searches together exclude a sneutrino LSP [15], at least if the majority of the CDM is the 
LSP. Thus we require in our analysis that the lightest neutralino be the LSP. 

Nevertheless there are new coannihilation processes to be considered when one or more 
of these 'wannabe' LSPs is almost degenerate with the hghtest neutralino. These include 
X - n, X - - ItL, X - - J^R, X - ^e,^i, X- X' - X^ coannihilations and all possible 
combinations ^. However, not all of these combinations are important as they are significant 
only in very small regions for a particular set of parameters. For this reason we do not 
include, for example, the sneutrino-slepton coannihilation in our calculations. 

We include in our subsequent calculations neutralino-slepton x — ^ [17, 7, 18], X ~ x' ~ 

^We neglect left-right (L-R) mixing for the first two generations of sfermions, so the right-handed selectron 
and smuon arc degenerate. Here and elsewhere, by 'right-handed' sfermion we mean the superpartner of 
right-handed fermion. 

squark LSP is possible only if \Ao\ is large, a possibility we do not study in this paper. 
''Again, because wc set Aq = here, squark coannihilations are not important, but see [16] for a calculation 
of neutralino-stop coannihilation. 
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[19], X ~ ^e,//, x' — and ~ ^ coannihilations^. The x' (co) annihilation rates can 
be derived from the corresponding x (co) annihilations by appropriate mass and coupling 
replacements. Details of our calculations are given in the Appendices. Following a summary 
of the relevant couplings in Appendix A, in Appendix B we update the neutralino-slepton 
coannihilation calculation of [17] to include L-R mixing. These are not very important at 
relatively low tan (5 but are potentially important for large values of tan (5. These improved 
coannihilation calculations were in fact already used in [21], but no details were given there. 
Appendix C provides chargino-slepton coannihilation processes, whilst Appendix D deals 
with neutralino-sneutrino coannihilation processes. 

4 Summary of Constraints and Review of the CMSSM 
Parameter Space 

We impose in our analysis the constraints on the MSSM parameter space that are provided 
by direct sparticle searches at LEP, including that on the lightest chargino x^'- "^x^ ^ 103.5 
GeV [22], and that on the selectron e: mg ^ 99 GeV [23]. Another important constraint is 
provided by the LEP lower hmit on the Higgs mass: niH > 114.4 GeV [24] in the Standard 
Model^. The hghtest Higgs boson h in the general MSSM must obey a similar hmit, which 
may in principle be relaxed for larger tan/3. However, as we discussed in our previous 
analysis of the NUHM [2], the relaxation in the LEP limit is not relevant in the regions of 
MSSM parameter space of interest to us. We recall that ruh is sensitive to sparticle masses, 
particularly mf, via loop corrections [25, 26], implying that the LEP Higgs limit constrains 
the MSSM parameters. 

We also impose the constraint imposed by measurements of 6 — > 57 [27], BR(S — > 
Xs-t) = (3.54 ± 0.41 ± 0.26) x 10"^, which agree with the Standard Model calculation 
BR{B X,-f)sM = (3.60 ± 0.30) x 10"^ [28]. We recall that the 6 ^ 57 constraint is more 
important for < 0, but it is also relevant for > 0, particularly when tan (5 is large, as we 
see again in this paper. 

We also take into account the latest value of the anomalous magnetic moment of the 
muon reported [29] by the BNL E821 experiment. The world average of = ^{g^ — 2) 
now deviates by (33.9 ± 11.2) x 10~^° from the Standard Model calculation of [30] using 
e+e~ data, and by (17 ± 11) x 10~^° from the Standard Model calculation of [30] based on 
^See [20] for recent work which includes all coannihilation channels. 

^In view of the theoretical uncertainty in calculating rUh, we apply this bound with three significant 
digits, i.e., our figures use the constraint mh > 114 GeV. 
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T decay data. Other recent analyses of the e+e~ data yield similar results. On some of 
the subsequent plots, we display the formal 2-a range 11.5 x 10~^° < < 56.3 x 10~^°. 
However, in view of the chequered theoretical history of the Standard Model calculations of 
a^, we do not impose this as an absolute constraint on the supersymmetric parameter space. 

As a standard of comparison for our NUHM analysis, we first consider the impacts of the 
above constraints on the parameter space of the CMSSM, in which the soft supersymmetry- 
breaking Higgs scalar masses are assumed to be universal at the input scale. In this 
mentioned in the Introduction, one may use the parameters (mi/2, mo, Aq, tan /3) and the sign 
of We assume for simplicity that Aq — 0, and plot in Fig. 1 the (mi/2, mo) planes for certain 
choices of tan /3 and the sign of ji. These plots are similar to those pubhshed previously [31], 
but differ in using the latest version of FeynHiggs [26] and the latest information on 
discussed above. 

The shadings and hues in Fig. 1 are as follows. The dark (brick red) shaded regions have 
a charged LSP, i.e. fi, so these regions are excluded. The 6 — > 57 exclusion is presented 
by the medium (dark green) shaded regions. The hght (turquoise) shaded areas are the 
cosmologically preferred regions with 0.1 < Q,^h^ < 0.3. The regions allowed by the E821 
measurement of at the 2-a level, 11.5 x 10"^*^ < < 56.3 x 10~^°, are shaded (pink) 
and bounded by sohd black hues. Only panel (a) and (d) have regions allowed by a^. The 
near- vertical (red) dot-dashed hues are the contours nih = 114 GeV, and the near- vertical 
(black) dashed line in panel (a) is the contour m^± = 103.5 GeV (though we do not plot 
this constraint in panels (b,c,d), the position of the chargino contour would be very similar). 
Regions on the left of these hues are excluded. 

We see in panel (a) of Fig. 1 for tan /3 = 10 and // > that all the experimental constraints 
are compatible with the CMSSM for mi/2 ~ 300 to 400 GeV and mo ~ 100 GeV, with larger 
values of mi/2 also being allowed if one relaxes the condition. In the case of tan/3 = 10 
and // < shown in panel (b) of Fig. 1, valid only if one discards the condition, the ruh and 
h ^ constraints both require mi/2 ^ 400 GeV and mo > 100 GeV. In the case tan /3 = 35 
and // < shown in panel (c) of Fig. 1, the 6 — > 57 constraint is much stronger than the ruh 
constraint, and imposes mi/2 ^ 700 GeV, with the Q.^h? constraint then allowing bands of 
parameter space emanating from mo ~ 600, 300 GeV. Finally, in panel (d) for tan /3 = 50 
and // > 0, we see again that the ruh and b ^ constraints are almost equally important, 
imposing mi/2 ^ 300 GeV for mo ~ 400 GeV. As in the case of panel (a), there is again a 
region compatible with the constraint, extending in this case as far as mi/2 ~ 800 GeV 
and mo ~ 500 GeV. 
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Figure 1: The CMSSM (mi/2,'n^o) planes for (a) tan/? = 10 and fi > 0, (b) tan/? = 10 
and < 0, (c) tan/5 = 35 and ^ < and (d) tan/3 = 50 and fi > 0, assuming Aq = 
0,Tnt = 175 GeV and mb{mb)g^]^ = 4.25 GeV. The near-vertical (red) dot-dashed lines are 
the contours ruh = 114 GeV as calculated using FeynHiggs [26], and the near-vertical (black) 
dashed line in panel (a) is the contour m^± = 103.5 GeV. The medium (dark green) shaded 
regions are excluded by b 57, and the light (turquoise) shaded areas are the cosmologically 
preferred regions with 0.1 < ^^^"^ — -^'^ dark (brick red) shaded regions, the LSP is 
the charged fi, so these regions are excluded. In panels (a) and (d), the regions allowed by 
the E821 measurement of at the 2-a level, as discussed in the text, are shaded (pink) and 
bounded by solid black lines, with dashed lines indicating the 1-a ranges. 
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5 Exploration of the NUHM Parameter Space 



Following our discussion of the CMSSM parameter space in the previous Section, we now 
discuss how that analysis changes in the NUHM. We extend our previous analysis [2] in two 
ways: (i) fixing tan/3 = 10 and > 0, but choosing different values of ^ and rriA, rather 
than assuming the CMSSM values, and (ii) varying tan P for representative fixed values of /i 
and niA- We make such selections for three projections of the NUHM, onto the (mi/2,mo) 
plane, the {/i, mA) plane and the (/x, M2) plane. 

5.1 The (mi/2,?7io) Plane 

Panel (a) of Fig. 2 shows the (mi/2, t^q) plane for tan /3 = 10 and the particular choice ji — 
400 GeV and = 400 GeV, assuming Aq = 0, m* = 175 GeV and mi,{rrn,)§M = 4.25 GeV as 
usual. Again as usual, the light (turquoise) shaded area is the cosmologically preferred region 
with 0.1 < Q^/i^ < 0.3. There is a bulk region satisfying this preference at mi/2 ~ 50 GeV to 
350 GeV, mo ~ 50 GeV to 150 GeV. The dark (red) shaded regions are excluded because a 
charged sparticle is fighter than the neutrafino. As in the CMSSM shown in Fig. 1, the fi is 
the LSP in the bigger area at larger mi/2, and there are fight (turquoise) shaded strips close 
to these forbidden regions where coannihilation suppresses the relic density sufficiently to 
be cosmologically acceptable. Further away from these regions, the relic density is generally 
too high. However, for larger mi/2 there is another suppression, discussed below, which 
makes the relic density too low. At small mi/2 and mo the left handed sleptons, and also the 
sneutrinos, become lighter than the neutrafino. The darker (dark blue) shaded area is where 
a sneutrino is the LSP. Within these excluded regions there are also areas with tachyonic 
sparticles. 

The near-vertical dark (black) dashed and light (red) dot-dashed fines in Fig. 2 are the 
LEP exclusion contours m^± > 104 GeV and ruh > 114 GeV respectively. As in the CMSSM 
case, they exclude low values of mi/2, and hence rule out rapid relic annihilation via direct- 
channel h and poles. The solid lines curved around small values of mi/2 and mo bound 
the light (pink) shaded region favoured by and recent analyses of the e^e~ data. 

A striking feature in Fig. 2(a) when mi/2 ~ 500 GeV is a strip with low Q,-xh^, which 
has bands with acceptable relic density on either side. The low-Q.^h?' strip is due to rapid 
annihilation via the direct-channel A,H poles which occur when m^^ = fnA/'^ — 200 GeV, 
indicated by the near-vertical solid (blue) line. Analogous rapid- annihilation strips have 
been noticed previously in the CMSSM [32, 21], but at larger tan/3 as seen in Fig. 1. There, 
they are diagonal in the (mi/2,?7^'o) plane, reflecting a CMSSM link between mo and m^ 
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Figure 2: Projections of the NUHM model on the (mi/2,mo) planes for tan (3 = 10 and 
(a) /i = 400 GeV and = 400 GeV, (h) fi = 400 GeV and = 700 GeV, (c) 
fi = 700 GeV and = 400 GeV and (d) fi = 700 GeV and ttla = 700 GeV, assum- 
ing Aq = 0,mt = 175 GeV and mb(mb)^j^ = 4.25 GeV. The near-vertical (red) dot-dashed 
lines are the contours nih = 114 GeV as calculated using FeynHiggs [26], and the near- 
vertical (black) dashed lines are the contours m^± = 103.5 GeV. The dark (black) dot-dashed 
lines indicate the GUT stability constraint. There are two such lines for each panel and only 
the areas in between are allowed by this constraint. The light (turquoise) shaded areas are the 
cosmologically preferred regions with 0.1 < Q^h"^ < 0.3. The dark (brick red) shaded regions 
is excluded because a charged particle is lighter than the neutralino, and the darker (dark 
blue) shaded regions is excluded because the LSP is a sneutrino. In panel (c) there is a very 
small medium (green) shaded region excluded by b ^ 57, at small mi/2. The regions allowed 
by the E821 measurement of at the 2-a level, as discussed in the text, are shaded (pink) 
and bounded by solid black lines. 
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that is absent in our implementation of the NUHM. The right-hand band in Fig. 2(a) with 
acceptable Q^/i^ is broadened because the neutralino acquires significant Higgsino content, 
and the relic density is suppressed by the increased W'^W~ production. Hereafter, we 
will call this the 'transition' band, which in this case is incidentally coincident with the 
right-hand rapid annihilation band ^. As mi/2 increases, the neutralino becomes almost 
degenerate with the second lightest neutralino and the lighter chargino, and the X — x' ~ 
coannihilation processes eventually push Q^h'^ < 0.1 when mi/2 ^ 700 GeV. We note that 
chargino-slepton coannihilation processes become important at the junction between the 
vertical bands in Fig. 2(a) and the neutralino-slepton coannihilation strip that parallels the 
— m~ boundary of the forbidden (red) charged-LSP region. 

There are two dark (black) dash-dotted lines in Fig. 2(a) that indicate where scalar 
squared masses become negative at the input GUT scale for one of the Higgs multiplets, 
specifically when either (mi(Mx)^ + n{MxY) < or {m2{MxY + //(Mx)^) < 0. One of 
these GUT stability fines is near-vertical at mi/2 ~ 600 GeV, and the other is a curved fine 
at mi/2 ~ 150 GeV, mo ~ 200 GeV. We take the point of view that regions outside either of 
these lines are excluded, because the preferred electroweak vacuum should be energetically 
favoured and not bypassed early in the evolution of the Universe, but a different point of 
view is argued in [13]. 

Thus, combining all the constraints, the allowed regions are those between the rrih line 
at mi/2 ~ 300 GeV and the stability line at mi/2 ~ 600 GeV, which include two rapid- 
annihilation bands, some of the transition band and the junction between the bulk and 
coannihilation regions around mi/2 ~ 350 GeV, mo ~ 150 GeV. If one incorporates also the 
putative a^j, constraint, only the latter region survives. We note however, that if the r data 
were used in the g — 2 analysis, the constraint from only excludes the lower left corner of 
the plane and large values of mi/2 and mo survive at the 2(7 level. 

Panel (b) of Fig. 2 is for // = 400 GeV and m^ = 700 GeV. We notice immediately that 

the heavy Higgs pole and the right-hand boundary of the GUT stability region move out 

to larger mi/2 ~ 850, 1050 GeV, respectively, as one would expect for larger m^- At this 

value of m^, the transition strip and and the rapid annihilation ('funnel') strip are separate. 

However the latter would be to the right of the transition strip and hence the Q.^h? bands on 

both sides of the rapid- annihilation strip that was prominent in panel (a) have disappeared, 

^As the neutralino acquires more Higgsino content, annihilation to W^W~ production increases, whilst 
fermion-pair production decreases (except for tt). Around mi/2 ~ 625 GeV, there is a threshold for hA, hH 
production, which decreases ^^h? to ~ 0.08, not far below the preferred range. However, the decrease in 
fermion production quickly raises Vt-^h^ again for larger mi/2- There is a very narrow stripe with fl^h^ < 0.1 
of width 5mi/2 ~ 2 GeV, which is not shown in the figure due to problems of resolution. 
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due to enhanced chargino-neutralino coannihilation effects. Panel (b) fias tfie interesting 
feature that there is a region of mi/2 ~ 300 to 400 GeV where mo = is allowed. As 
discussed in [33] this possibility, which would be favoured in some specific no-scale models of 
supersymmetry breaking, is disallowed in the CMSSM. The small-mo region is even favoured 
in this variant of the NUHM by the putative constraint. 

Panel (c) of Fig. 2 is for fj, — 700 GeV and m^ = 400 GeV. In this case, we see that 
the rapid- annihilation strip is back to mi/2 ~ 500 GeV, reflecting the smaller value of m^, 
whereas the transition band has separated off to large mi/2, reflecting the larger value of 
fjL. However, this band is excluded in this case by the GUT stability requirement. GUT 
stability also excludes the possibility that mo = 0. In this case, the putative constraint 
would restrict one to around mi/2 ~ 400 GeV and mo ~ 100 GeV. 

Finally, panel (d) of Fig. 2 is for // = 700 GeV and m^ = 700 GeV. In this case, the 
rapid-annihilation strip has again moved to larger mi/2, related to the larger value of uia, 
and the transition band at large mi/2 is again excluded by the GUT stability requirement. 
The bulk region has disappeared in this panel, reflecting the fact that the values of /i and m^ 
here has strayed away from their values in the bulk region for the CMSSM. GUT stabihty no 
longer excludes mo = 0, and this possibility would be selected by the putative constraint. 

We now turn, in Fig. 3, to the impact of varying tan/3, keeping // = 400 GeV and 
m^ = 700 GeV. For convenience, panel (a) reproduces Fig. 2(b) with tan/3 = 10. In general 
as tan/3 is increased, two major trends are visible. One is for the region excluded by the 
requirement that the LSP be neutral to spread up to larger values of mo, and the other is 
for the constraint to move out to larger values of mi/2 and mo. 

Specifically, we see in panel (b) of Fig. 3 for tan P — 20 that, whereas the heavy Higgs 
pole at mi/2 ~ 850 GeV essentially does not move, the fi LSP and coannihilation strip lying 
above the excluded charged-LSP region rise to larger mo. This has the effect of excluding 
the mo = option that was present in panel (a). At low mi/2, the nih constraint is stronger 
than the GUT stability and other constraints. The constraint would allow a larger range 
of mi/2 than in panel (a), extending up to ~ 550 GeV. 

Continuing in panel (c) of Fig. 3 to tan /3 = 35, we see that the minimum value of mo has 
now risen to ~ 200 GeV. We also see that the 6 — > 57 constraint is now important, enforcing 
mi/2 > 300 GeV in the region preferred by the relic density. Because of this and the irih 
constraint, the GUT stability constraint is now irrelevant at low mi/2, whereas at high mi/2 
it has vanished off the screen, and is in any case also irrelevant because of chargino-neutralino 
coannihilation. The constraint would now allow part of the cosmological band on the left 
side of the rapid- annihilation strip. These trends are strengthened in panel (d) of Fig. 3 for 
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Figure 3: The NUHM {mi/2,mQ) planes for (a) tan/5 = 10 , (b) tan/3 = 20, (c) tan/5 = 35 
and (d) tan/5 = 50, for n = 400 GeV, = 700 GeV, assuming Aq = 0, = 175 GeV and 
mb^mbjs^M = 4.25 GeV. The shadings and line styles are the same as in Fig. 2. 
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tan/3 = 50, where we see that mo,mi/2 ^ 400 GeV because of the 6 — > 57 constraint, and 
would allow mo ^ 1100 GeV along the transition band. 

5.2 The (//,m^) Plane 

We now analyze the range of possibilities in the (/i, m^) plane for various fixed choices of 
mi/2 and mo, first choosing tan/? = 10. Panel (a) of Fig. 4 displays the (/i, m^) plane for 
^1/2 = 300 GeV, mo = 100 GeV. As we saw in Fig. 1(a), there is a CMSSM point with 
/i > that is compatible with all the constraints for these values of mi/2 and mo. The 
corresponding CMSSM point for |U < in Fig. 1(b) is, however, incompatible with the m^ 
and 6 — 57 constraints, as well as the putative constraint. The CMSSM equivalent points 
are shown as crosses in panel (a) of Fig. 4, and have (/i, itla) — (±390, 450) GcV. 

As usual, there are dark (red) regions where there is one or more charged sparticle 
lighter than the neutralino x so that x is no longer the LSP. First, there are 'shark's teeth' 
at ~ 300 GeV, m^ ^ 300 GcV in panel (a) of Fig. 4 where the fi is the LSP. At small 
particularly at small m.A when the mass difference rn^ — m^ is small, the mass is 
driven small, making the fi the LSP At even smaller |^|, however, the lightest ncTitralino 
gets lighter again, since m^ ~ /i when fi < Mi ~ 0.4mi/2 In addition, there are 
extended dark (red) shaded regions at large \fi\ where left-handed sleptons become lighter 
than the neutralino. However, the electron sneutrino and the muon sneutrino (which 
are degenerate within the approximations used here) have become joint LSPs at a slightly 
smaller \n\. Since the possibility of sneutrino dark matter has been essentially excluded by 
a combination of 'u counting' at LEP, which excludes ma < 43 GeV [34], and searches for 
cold dark matter, which exclude heavier i> weighing < 1 TeV [15], we still demand that the 
LSP be a neutralino x- The darker (dark blue) shaded regions are where the sneutrinos are 
the LSPs, and therefore excluded. 

To explain the possible spectra better, we plot in Fig. 5 some sparticle masses as functions 
of II, with other parameters fixed. In particular, we plot the neutralino mass m^ (dark solid 
curve), the chargino mass m^± (dark dashed curve), the lighter stau mass m~ (hght sold 
curve), the right-handed selectron mass m~^ (light dashed curve), and the sneutrino masses 
m~ and m~ (thin solid and dashed curves respectively). We have omitted the curve for 
m~^, which is very similar to the sneutrino curves, for reasons of clarity. In both panels of 

®Note that the cr and JIr masses are also driven to small values, along with the tr, and in fact there are 
small regions where the degenerate cr and Jir become the LSP. 

'^'^We note that regions in Figs. 4(a,b) and 7(b) at small and niA are excluded by the LEP slectron 
search [23]. However, these regions are all excluded by other limits, such as GUT stability, the LEP chargino 
limit [22] and b — > S7, and so are not shown separately. 
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Figure 4: The NUHM (/i, m^) planes for tan (3 = 10, (a) mo = 100 GeV and mi/2 = 
300 GeV, (b)mo = 100 GeV and mi/2 = 500 GeV, (c) uiq = 300 GeV and mi/2 = 300 GeV 
and (d) uiq = 300 GeV and mi/2 = 500 GeV, assuming Aq = 0, = 175 GeV and 
^bi'^bJs^M — 4-25 GeV. The shadings and line styles are the same as in Fig. 2. 
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Fig. 5, we have fixed most parameters as in Fig. 4(a). Fig. 5(a) is for m.A — 200 GeV, and 
we see that the hghter stau and the sneutrinos are hghter than the neutrahno for very large 
values of As |//| is decreased, there is a small window in which i/g is the LSP. The curious 
behaviour of the stau mass curve at large |//| is due to a large and negative value for S which 
drives down while driving up (cf. Eq. (7)). For large values of the lighter stau 
is mostly left-handed, whilst it is mostly right-handed for smaller |//| < 1000 GeV. When 
|//| ^ 400 GeV, both the stau and the right-handed selectron are hghter than the neutrahno 
until is very small and the neutrahno becomes Higgsino-like, with its mass scaling as /i. 
The 'shark's teeth' in Fig. 4(a) correspond to the range around |//| ~ 200 GeV where the 
lighter stau is the LSP, with the right-handed selectron shghtly heavier. 

In Fig. 5(b), we have fixed = 600 GeV. Here, one sees a similar pattern of masses 
at large |//|. However, at small the lighter stau and right-handed selectron are much 
heavier, so the neutrahno remains the LSP. This reflects the fact that uia — 600 GeV is 
above the tips of the 'shark's teeth' in Fig. 4(a). It is clear that the sizes of these 'shark's 
teeth' must depend sensitively on the NUHM parameters, as seen in the other panels of 
Fig. 4 and subsequent flgures. 

Returning to panel (a) of Fig. 4, we see strips adjacent to the Ve,n LSP regions, where 
neutralino-sneutrino coannihilation is important in suppressing the relic density to an accept- 
able level. Appendix D presents details of our calculations oi i> — x coannihilation channels. 
We consider here only z/g and U/j,, for which the effects of mg, and L-R mixing in the slep- 
ton mass matrix are negligible, leading to some simpliflcations compared to the Ur (whose 
inclusion would have little effect on our flgures) . The inclusion of neutralino-sneutrino coan- 
nihilation in panel (a) squeezes inwards slightly the coannihilation strips. As we see later, 
the effect of neutralino-sneutrino coannihilation is more noticeable at larger tan/3. Again as 
usual, the hght (turquoise) shaded regions are those for which 0.1 < fl^h'^ < 0.3. 

The thick cosmological region at smaller /i in panel (a) corresponds to the 'bulk' region 
familiar from CMSSM studies. Note that this region is squeezed outward near the 'shark's 
teeth' by the neutralino-slepton coannihilation. Extending upward in ttia from this bulk 
region, there is another light (turquoise) shaded band at smaller |//|. This is the transition 
band, where the neutrahno gets more Higgsino-content and the annihilation to W'^W~ 
becomes important, yielding a relic density in the allowed range, as happens in the focus- 
point region [35] in the CMSSM. For smaller |//|, the relic density becomes too small due to 
X ~ x' ~ coannihilations, and the chargino-slepton annihilations described in Appendix 
C must be taken into account where this strip meets the neutralino-slepton coannihilation 
strip discussed earher. The LEP hmit on the chargino mass excludes a strip at even smaller 
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Figure 5: The masses (dark solid), 'rn^± (dark dashed), m~ (light solid), (light 
dashed), (thin solid) and m~ (thin dashed) as functions of fi for tan (3 = 10, mi/2 = 
300 GeV, mQ = 100 GeV for Jafm,^ = 200 GeV and (b) = 600 GeV, assuming Aq = 
0, mt = 175 GeV and nibimb)^ = 4.25 GeV 
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where the rehc density would again have come into the cosmologically preferred region. 
There are also horizontal bands of acceptable rehc density when uia ~ 250 GeV, that are 
separated by strips of low relic density, due to rapid annihilation through the A (indicated 
by sohd (blue) lines) and H poles. 

Underlying the cosmological regions are dark (green) shaded regions excluded by 6 — > 57, 
which are more important for // < 0, as expected from previous analyses. Also important 
is the nih constraint, calculated using the FeynHiggs programme [26], which excludes the 
option < in panel (a) of Fig. 4. The putative constraint would also exclude the < 
half-plane, while allowing all of the // > parameter space for this particular choice of tan/5, 
mi/2 and mo. 

The darker (black) dot-dashed hues in Fig. 4(a) indicates where one or the other of the 
Higgs mass-squared becomes negative at the input GUT scale. We see that these constraints 
exclude much of the cosmological region still permitted, apart from the 'bulk' region and 
part of the 'transition' region for small // > 0. 

We conclude from Fig. 4(a) that moderate positive values of < 700 GeV are favoured 
and niA is unlikely to be very small, though there is a very small allowed region below 
the rapid A, H annihilation strip where niA ~ 230 GeV. The cosmological relic density hes 
within the range favoured by astrophysics and cosmology in a large fraction of the remaining 
NUHM parameter space for > 0, generahzing the CMSSM point that is indicated by the 
cross. 

The notations used for the constraints illustrated in the other panels of Fig. 4 are the 
same, but the constraints interplay in different ways. In panel (b), we have chosen a larger 
value of mi/2. In this case, the dark (red) shaded 'shark's teeth' at moderate |//| and small 
niA where the LSP is charged have expanded greatly, and are flanked by pale (turquoise) 
shaded regions where neutralino-slepton coannihilation produces an acceptable relic density. 
On the other hand, the dark (red) shaded and the darker (dark blue) shaded regions at 
large |//| have moved out of the panel displayed, and one sees only parts of the adjacent 
coannihilation strips. The relatively large value of mi/2 keeps the rate of 6 — > 57 under 
control unless is small and /i < 0. The chargino constraint is similar to that in panel 
(a), whereas the nih constraint is irrelevant due to the large value of mi/2- The putative 
constraint would allow only a very small region in this panel, but without any cosmological 
preferred region. Finally, we observe that the GUT stability constraint now allows larger 
values of \n\ < 1000 GeV and tua ^ 300 GeV. 

We now turn to panel (c) of Fig. 4, which is for tan/3 = 10, mi/2 = 300 GeV and 
mo — 300 GeV. In this case, the 'shark's teeth' have disappeared, as have the neutrahno- 
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slepton strips at large |//| (due to the large value of mo). Negative values of fj, are excluded 
partially by mh and by 6 ^ 57, as well as by the putative constraint, which would permit 
a strip with // > 0. GUT stabihty enforces // ^ 800 GeV, and also provides a lower limit 
on niA that is irrelevant because of the other constraints. In contrast, in panel (d) of Fig. 4 
for mi/2 = 500 GeV and mo = 300 GeV, we see a similar 'cruciform' pattern for the regions 
allowed by cosmology, but 6 — > 57 has only rather limited impact for // < and mh is 
irrelevant: there is no region consistent with the putative constraint. In this case, m^ 
could be as small as the ~ 300 GeV allowed by the GUT stabihty constraint, and could 
be as large as ~ 1100 GeV. 

We now discuss the variation with tan (3 shown in Fig. 6, keeping fixed mi/2 = 300 GeV 
and mo = 100 GeV, and starting (for convenience) in panel (a) with the case tan P — 10 
shown previously in panel (a) of Fig. 4. Increasing tan /3 to 20 in panel (b), we note that the 
'shark's teeth' are somewhat expanded, whereas the regions at large and/or ttia where the 
hghtest neutralino is no longer the LSP change in shape. The allowed cosmological region 
is bounded at larger m^ by one where the LSP is the fi, and at larger |//| by one where 
the LSP is a sneutrino: these different regions had shapes similar to each other in panel (a). 
As usual, neutralino-slepton coannihilations - as calculated in Appendix B - and neutralino- 
sneutrino coannihilations - as calculated in Appendix D - suppress the relic density close to 
these boundaries, and the relic density is within the preferred range over most of the region 
outside the 'shark's teeth' and the chargino exclusion for m^ > 2m^. As in previous cases, 
6 — > 57 and m^ together exclude the option /i < 0, which would also be disfavoured by 
a^. The fi > region surviving all the constraints when tan/3 = 20 hes between the GUT 
stabihty constraint at ~ 600 GeV and 6 — > 57 constraint around // ~ 350 GeV, and has 
130 GeV^ ruA ^900 GeV. 

When tan/3 = 35, in panel (c) of Fig. 6, we find no consistent electroweak vacuum for 
a large (polka-dotted) region with // < 0, and the condition that the LSP not be the fi 
excludes a broad swathe with small // > 0. The condition that the LSP not be a sneutrino 
provides an upper bound ^ 1500 GeV. The nih constraint imposes m^ ^ 120 GeV and 
the b ^ constraint is irrelevant. The GUT stability constraint allows only a narrow sliver 
of ~ 700 GeV for 120 GeV < mA < 200 GeV, where the lower bound comes from the LEP 
Higgs search. This reach is not compatible with a^. Here, the predicted discrepancy with 
the Standard Model is in excess of that allowed at the 2a level by the g — 2 experiment. 

In the case of tan /3 = 50, shown in panel (d) of Fig. 6, the problem with the non-existence 
of a consistent electroweak vacuum extends to most of the /i < half-plane, as shown by 
the polka dots. The regions excluded because the LSP is either the fi or a sneutrino hem in 
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Figure 6: T/ie NUHM (/i, m^) planes for uiq = 100 GeF anc? mi/2 = 300 GeV, for (a) 
tan/? = 10, (b) tan/5 = 20, (c) tan/5 = 35 and (d) tan/5 = 50, assuming Aq = 0,mt = 
175 GeV and mb{mb)g^M = 4.25 GeK T/ie shadings and line styles are the same as in Fig. 2, 
and there is no consistent electroweak vacuum in the polka-dotted region. 
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a small region with // > 0. The region that survives all the experimental and cosmological 
constraints is similar to that for tan (3 — 35, but extends to somewhat larger ttia and broader 
II. Once again, these regions are not compatible with the g — 2 result. 

Another series of plots for different values of tan /3 is shown in Fig. 7, this time for the 
fixed values mi/2 = 500 GeV and mo = 300 GeV. All of these have the 'cruciform shape' of 
cosmological region famihar from panels (c, d) of Fig. 4. The constraint is irrelevant for 
this larger value of mi/2, but much of the parameter space for // < is excluded by & — > 57, 
particularly for larger tan/3. Likewise, a^j, consistently favours > 0. One of the GUT 
stabihty constraints requires |//| ^ 1200 GeV, the exact value increasing slightly with tan/3. 
The other GUT stability constraint requires niA <^ 300 GeV for tan /3 = 10, but weakens for 
larger tan /3, becoming irrelevant when tan /3 > 35. As in Fig. 6, the absence of a consistent 
electroweak vacuum also becomes a problem for // < when tan/3 > 35. Finally, we see in 
panel (d) of Fig. 7 that small values of positive fi are disallowed because the LSP becomes 
the fi. As in Fig. 6, values of niA as low as allowed by the Higgs search are consistent with 
all the constraints for tan/3 > 35, as long as fj, has a suitable positive value. Unhke the 
case shown, in Fig. 6, here the entire /i > half-plane is within the 2a range for a^^ when 
tan P — 35 and 50. 

5.3 The (//, M2) Plane 

Panel (a) of Fig. 8 displays the (/i, M2) plane for the choices tan/3 = 10, mo = 100 GeV 
and niA = 500 GeV. We restrict our attention to the region allowed by the GUT stability 
constraints, which is roughly triangular, extending from the origin to vertices at (/x, M2) = 
(±1500, 600) GeV. Within this region, there are two large triangular regions with M2 ^ 
300 GeV and either sign of /j, that are excluded because the LSP is not the lightest neutralino. 
Also, a band around the = axis is excluded by the LEP chargino constraint, regions with 
M2 ^ 250(300) GeV for fj, > (<)0 are excluded by the LEP Higgs constraint, and most of the 
surviving < region is eroded by the 6 — > 57 constraint. In the lower corners (for both signs 
of /i) the LSP is a sneutrino and we see the effects of x — coannihilation running alongside 
these regions. However, these regions are in conflict with the GUT stability constraint. 

Much of the remaining area of the plane is consistent with the cosmological relic den- 
sity constraint, mainly along strips where neutralino-slepton is important. These terminate 
around M2 ~ 450 GeV, where rapid annihilations through the direct-channel A, H reso- 
nances cut down the relic density. The putative constraint would favour /i > and 
M2 < 350 GeV. 
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Figure 7: T/ie NUHM (/i, m^) planes for uiq = 300 GeF anc? mi/2 = 500 GeV, for (a) 
tan/? = 10, (b) taii/S = 20, (c) tan/3 = 35 and (d) tan/5 = 50, assuming Aq = 0,mt = 
175 GeV and mb{mh)g^j^f = 4.25 GeV. The remaining shadings and line styles of Fig. 2 are 
used here, and there is no consistent electroweak vacuum in the polka-dotted region. 
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Figure 8: The NUHM [fi, Ma) planes for tan (3 = 10, (a) mo = 100 GeV and = 500 GeV, 
(h) mo = 100 GeV and = 700 GeV, (c) mo = 300 GeV and m^ = 500 GeV and_(d) 
mo = 300 GeV and m^ = 700 GeV, assuming Ao = 0,mt = 175 GeV and mb{mb)^M = 
4.25 GeV. The shadings and line styles are the same as in Fig. 7. 
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The above constraints interplay analogously in panel (b) of Fig. 8, for mo = 100 GeV and 
niA — 700 GeV. We note, however, that the GUT stabihty hmit has risen with niA to above 
800 GeV, and that the tips of the non-neutralino LSP triangles have also moved up shghtly. 
The neutralino-slepton coannihilation strip now extends up to M2 ~ 640 GeV, where it is 
cut off by rapid A, H annihilations. We also see the 'transition' band at lower M2. In fact, 
this region with // > is the one favoured by all constraints including a^. 

In panel (c) of Fig. 8, for mo = 300 GeV and m^ = 500 GeV, the triangular GUT 
stability region is very similar to that in panel (a), whereas the non-neutralino LSP triangles 
have moved to significantly higher M2, reflecting the larger value of mo. The regions where 
the cosmological relic density falls within the preferred range are now narrow strips in the 
neutralino-slepton coannihilation regions, on either side of the rapid A, H annihilation strips, 
and the 'transition' bands. This tendency towards 'skinnier' cosmological regions is also 
apparent in panel (d), where mo = 300 GeV and vriA — 700 GeV are assumed. In both 
these panels, the b ^ sj constraints disfavours /i < and small m^, whilst the putative 
constraint would favour /i > and small m^- 

We now explore the variation of these results with tan /3, as displayed in Fig. 9 for the 
choices mo = 100 Gev and uia — 500 GeV. Panel (a) reproduces for convenience the case 
tan P — 10 that was shown also in panel (a) of Fig. 8. When tan /3 = 20, as seen in panel (b) 
of Fig. 9, we first notice that the GUT stability region now extends to larger M2. Next, we 
observe that the triangular non-neutralino LSP regions have extended down to lower M2. As 
one would expect on general grounds and from earlier plots, the mh constraint at low M2 is 
weaker than in panel (a), whereas the b ^ constraint is stronger, ruling out much of the 
otherwise allowed region with /i < 0. The allowed region with // > is generally compatible 
with the putative a^u constraint. 

Panel (c) for tan (3 — 35 again shows the feature that no consistent electroweak vacuum 
is found over much of the half-plane with // < 0. The GUT stabihty constraint now allows 
M2 ^ 1000 GeV. The non-neutrahno LSP region is no longer triangular in shape, but now 
requires // ^ 800 GeV. This happens as the regions with hght left-handed slepton at low 
M2 meet with the ones with light right-handed slepton at higher M2. There is a minuscule 
allowed region for < 0. The residual regions with relic density in the preferred range for 
/I > are limited to strips in the neutralino-slepton coannihilation regions, and on either 
side of the rapid A, H annihilation strip. All this preferred region is compatible with the 
putative constraint. 

A rather similar pattern is visible in panel (d) of Fig. 9 for tan/3 = 50. The most 
noticeable differences are that the absence of a consistent electroweak vacuum is even more 



25 



tan p = 10 , m„ = 100 GeV, = 500 GeV tan p = 20 , m„ = 100 GeV, = 500 GeV 




(GeV) (GeV) 

Figure 9: The NUHM (/x, M2) planes for mo = 100 GeV and = 500 GeV, for (a) tan (3 = 
10, (h) tan/3 = 20, (c) tan/3 = 35 and (d) tan/? = 50, assuming Aq = 0,mt = 175 GeV and 
mb{mb)s^M = 4.25 GeV. The shadings and line styles are the same as in Fig. 7. 
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marked for > 0, and that the putative constraint would suggest a lower limit M2 ^ 
350 GeV, whereas lower values would have been permitted in panel (c) for tan/3 = 35. 

6 Conclusions and Open Issues 

We have provided in this paper the tools needed for a detailed study of the NUHM, in the 
form of complete calculations of the most important coannihilation processes. However, in 
this paper we have only been able to scratch the surface of the phenomenology of the NUHM. 
Even this exploratory study has shown that many new features appear compared with the 
CMSSM, as results of the two additional parameters in the NUHM, but much more remains 
to be studied. For example, we have not studied the NUHM with a nonzero trilinear couphng 
Aq. Nevertheless, some interesting preliminary conclusions about the NUHM can be drawn, 
though many questions remain open. 

The lower hmit on the LSP mass in the CMSSM has been much discussed, and it is 
interesting to consider whether this could be greatly reduced in the NUHM. The top panel 
of Fig. 10 compiles the bounds on for the various sample parameter choices explored 
in this paper. We focus on the case /i > 0, which is favoured by rrih and 6 — > 57 as well 
as the dubious — 2 constraint. The solid (black) line connects the lower limits on 
for the particular choice ji — 400 GeV and = 700 GeV shown for the four choices of 
tan /3 in Fig. 3. This lower hmit is provided by ruh for tan /3 = 10, 20, but by 6 ^ 57 for 
tan 13 — 35, 50. The CMSSM lower bound on is indicated by a thick (blue) line. We see 
that this CMSSM lower bound on m-x^ is similar to that in the NUHM when tan/3 = 10, but 
weaker when tan /3 = 50 because of the different behaviour in the CMSSM of the 6 — > 57 
constraint in this case. 

For tan (3 — 10, we also indicate by different black symbols the lower bounds on for 
the other choices of {ji^mA) studied in Fig. 2, namely {ji^mA) — (400,400) GeV (square), 
{n,mA) = (700,400) GeV (diamond) and (/i,TOA) = (700,700) GeV (star). The two latter 
lower bounds are significantly higher than in the default case (//, m^) = (400, 700) GeV, 
due to the impacts of the GUT stability condition and the Q,y.h? constraint, respectively. 
In none of the NUHM examples studied was the lower limit on relaxed compared with 
the CMSSM, though this might be found possible in a more complete survey of the NUHM 
parameter space. 

We also show in the top panel of Fig. 10 the upper bounds on found in the NUHM for 
the various parameter choices explored in this paper. The (red) dot-dashed line is the upper 
bound for ji — 400 GeV and = 400 GeV, as obtained without imposing the putative 
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Figure 10: Top panel: Bounds on the LSP mass found in the NUHM for the particular 
parameter choices displayed in Figs. 2 and 3. The solid (black) line is the lower limit on 
for {fi,mA) = (400,700) GeV, The (red) dotted line is the upper bound for the same values 
of {fj,,mA), but without imposing the putative — 2 constraint, and the (green) dashed line 
shows how this upper bound would be strengthened if the g^ — 2 constraint were imposed. The 
symbols shown for tan /3 = 10 mark the lower ( g^~2 upper, overall upper) limits on for 
(/i, m^) = (400,400) GeV (squares), {^,mA) = (700,400) GeV (diamonds) and {fi,mA) = 
(700, 700) GeV (stars). The thick (blue) lines correspond to GMSSM upper and lower limits. 
Bottom panel: Lower bounds on rriA for {mi/2,fnQ) = (300, 100) GeV (solid black line) and 
(mi/2,'^o) = (500,300) GeV (dashed red line). The extra symbols for tan f3 = 10 are lower 
bounds for {mi/2,mQ) = (500,100) GeV (upper point) and (300,300) GeV (lower point). 
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— 2 constraint. The upper bound on is attained in the rapid-annihilation strip where 
< niA/^. Also shown in Fig. 10, as a (green) dashed line, is the strengthening of this 
upper bound that would be found if the — 2 constraint were imposed. This constraint 
would not strengthen the upper limit for tan /3 > 35. We also show as a thick (blue) line the 
CMSSM upper limit on m^^, implementing the — 2 constraint. We see that it is similar 
to that in the NUHM for tan /3 = 10, 20, but is weaker for higher tan p. The CMSSM upper 
bound would be much weaker still if the g^ — 2 constraint were relaxed, because of the 
different behaviours of the rapid-annihilation strips in the NUHM and CMSSM, as seen by 
comparing Figs. 1 and 3. 

For tan /3 = 10, we also show in the top panel of Fig. 10 the upper bounds that hold with 
and without g^ — 2 for the other three choices of (/x, m^) made in Fig. 2, using the same 
symbols as for the lower bounds (squares, diamonds and stars, respectively). The largest 
upper bound is for the choice (/x, m^) = (700, 700) GeV, for which there is an extension of 
the allowed region along the coannihilation tail, extending up to the GUT stabihty hne. 

In the second panel of Fig. 10, we summarize the lower limits on found in the 
NUHM models that we have studied. The solid (black) hne is for the choice {mi/2,mo) — 
(300, 100) GeV displayed previously in Fig. 7. The lower bounds shown for tan (3 — 10, 20 
are compatible with g^ — 2, but we find that no region would be allowed hy g^ — 2 in the 
cases tan /3 = 35, 50, as seen in panels (c) and (d) of Fig. 7. The (red) dashed line in the 
second panel of Fig. 10 is for the choice (mi/2,mo) = (500,300) GeV shown previously in 
Fig. 6. In this case, we find that no region would be allowed by g'^ — 2 for tan /3 = 10, whereas 
the lower bounds for tan /3 = 20, 35, 50 are compatible with g^ — 2. The extra symbols for 
tan /3 = 10 correspond to the other choices (mi/2, ?7io) = (500, 100) GeV and (300, 300) GeV 
shown in panels (b) and (c) of Fig. 4, at = 308 GeV and = 220 GeV, respectively. 
We note that neither of these choices satisfies the g — 2 constraint. 

An interesting feature of the second panel of Fig. 10 is the fact that the lower bounds coin- 
cide for tan /3 > 35, and correspond to the lower limit established by direct searches at LEP. 
For comparison, we note that in the CMSSM for (mi/2, mo) = (300, 100) GeV one would have 
ruA = 449, 424, 377, 315 GeV for tan (5 = 10, 20, 35, 50, whilst for (mi/2, mo) = (500, 300) GeV 
one would have m^ = 762,720,639,526 GeV for tan/3 = 10,20,35,50. We conclude that 
the NUHM allows m^ to be considerably smaller than in the CMSSM, particularly at large 
tan 13. 

This brief survey is no substitute for a detailed study of the NUHM. However, we have 
provided in the Appendices the technical tools required for such a study, and in this Section 
we have presented some preliminary observations based on a cursory exploration of the 
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NUHM. This has already provided some interesting indications, for example that it may 
prove difficult to relax significantly the CMSSM lower bound on the LSP, but that itia may 
be greatly reduced. In our view, it would be interesting to pursue these questions more 
deeply in a detailed study of the NUHM, a project that hes beyond the scope of this work. 
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Appendix A: Couplings used in the Calculations 

Here we list the couplings used in the calculations. For clarity, we have not written fac- 
tors such as i, 7's and momenta. These are taken into account in the calculations of the 
amplitudes squared in the other Appendices. 

Couplings for chargino-slepton coannihilation: 
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Couplings for neutralino-sneutrino coannihilation: 
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Ch-h-h = -3c/2mz/ (2 cos 6*^.) cos(2q;) sin(/3 + a) 
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^i/2((l/2)/cos2^^(-cos(2Q;))) 
gl/2{{l/2)/cos^e^ cos(2«)) 
^2V2sin(2a)((l/2)/cos2^^) 
difrizj (2 cos cos(2/3) cos(/3 + a) 
—g2fnzl (2 cos ^^k) cos(2/3) sin(/3 + a) 
52V2((l/2)/cos2^^^,cos(2/3)) 
gi sin(/9 — a)/{2 cos 6*^1.) 
— cos(/9 — a)/{2 cos ^v^.) 
(?2/(2cos^,^)(iVi4Ari4 - ATisATia) 

-^2/(2 cos ^^)(iVi47Vl4 - NMs) 

-g2/V2 

-92/V2 

-c/2/ V2(sm 6^2 tan ^i^iVn) 
-^2/V2(cos^e2 tan^^TV^i) 
-g2/^(cosee(-Ni2-ta.ne^Nii)) 
-^2/v^(- sin^e(-A^i2 - tan^^A^n)) 
^2(-l/y2iVi4V.2 + iVi2Ka) 
^2(l/V27Vi3[/,2 + 7Vi2t/a) 

g2/2((Ni3(Ni2 - Nil tan ^^^) + Nis(Ni2 - A^ii tan ^^^)) sin a 
+ {Nu{Ni2 - Nil tanew)Ni4{Nu - Nn tanO^)) cos a) 
^2/2((iVi3(iVi2 - A^ii tan 6^,) + 7Vi3(iVi2 - Na tan 6^)) sin a 
+{Ni4{Ni2 - Nu tan 9^) + Nu{Ni2 - Nn tan 0^)) cos a) 
92/'2((Ni3{Ni2 - Nil tan + A^i3(A^i2 - A^ii tan ^vk))(- cos a) 
+ (A^i4(iVi2 - Niitd.ne„)Ni4{Ni2 - NntanOw)) sin a) 
52/2((Ari3(A^i2 - A^iitan^^) + Ni-i{Ni2 - Nn tan^^))(-cosQ;) 
+{Ni4{Ni2 - Nil tan 9^) + Nu{Ni2 - Nn tan 0^^)) sin a) 
92/2{{Ni3{Ni2 - Nil tan 9^,) + Nis{Ni2 - Nn tan 9^)) sin /3 
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+{Nu{Ni2 - Nateine^)Ni4{Ni2 - A^n tan^;^))(- cos/3)) 
g2/2{{N,3{Ni2 - Nu tanO^) + Nis{N^2 - iV,i tan^^))(- sin/?) 
+ {Ni4{Ni2 - Nu tan Oy^) + Nu{Ni2 - Na tan^^^)) cos/3) 
-^2/v^sin(2/3) 

-92{NuVny/l/2{Ni2 + Nn tane^)Vi2) cos/3 
-g2{NMi - yi72(iVi2 + Nntaney,)Ui2) sin/3 

Appendix B: Neutralino-Slepton Coannihilation with L- 
R Mixing 

We list below some modifications and additions to the slepton-slcpton and ncutralino-slepton 
coannihilation channels previously given in the Appendix of [17]. There the lighter sleptons 
were assumed to be pure partners of the right-handed leptons, which is a very good approxi- 
mation for CMSSM when tan /3 is small. Here, we include L-R mixing and denote the hghter 
sleptons by Ii, denoting £ = e,iji,r. The obvious change needed is then replacing fi,, in [17] 
with ii. Furthermore, the presence of the left component entails some modifications in the 
couplings and opens up some new channels. We note that L-R mixing had already been 
included in some of our previous work [21, 2]. Our convention for the slepton mixing angle 
is 

( i^\ = ( ^^^^^ "^"^^ )( i^) 

[£2) [-sm0e COS0, J [ in J ^ ' 

The sign of m^o must be kept in all of the formulae below. 

£iil — ^ W+W- 

There are two additional channels: 

III. i-channel z?^ exchange 

IV. point interaction 

The couplings /i, /2 and /s are modified, while /e remains the same. 

/i = {— g2mw cos{P — a)){g2mz/ cos 9w{{i/2 — sin^ 9w) cos{P + a) cos^ Oi 
+ sin^ 9w cos(/3 -|- a) sin^ Oi) + g2Tnny{m^ cos /3)(— cos a) 
—g2'>Ti£/{mw COS (3) {—Ai COS a — //sin a) sin cos ^^) 
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'^xl-xT-H+ 

^X^-X--H+ 



/2 = {— g2mw sm(P — a))(g2mz/ cos9w((—'i-/'2 + sm'^ 6w) sm(P + a) cos^ 6i 
— sin^ 9w sin(/5 + a) sin^ 9e) + g2m'j / {rriw cos /3) sin a 
—g2mt/ {jUw cos /5) (^4^ sin a — cos a) sin 6*^ cos 9() 

h = (-^2Cose,/V2)2 

/4 = -gl cos^ 61^/2 

/s = (5-2 cos 6'i„.)(-5(2/ COS 6'iy(sin^ 6*^^ - cos^ 6*^/2)) 



rrir. — u 



TiyxTiy = {12mi - Amis + s^)/ {Ami: 
TixTui = (-6m|m^ - 20m|m^ - 6m^ + m|s + 2m|m^s + 5m^s + 4m|m^i 

+Am^t + Sm^^m^u + Sm^u — 2m^^su — 2m^su — Am^tu — 2m^u 

+su'^)/{2ml,{mjj - s)(m?^ - u)) 
TiixTui = {-Gmj^ml - 20m|mf,, - 6m^ + m|^s + 2m|m^s + 5m^s + Am\mlt 

-\-Am1t + Smj^mlu + Smlu — 2rnj^su — 2'mlsu — Aw^tu — 2rri^u'^ 

+su^)/{2ml,{ml - s)(m?^ - u)) 
TiuxTiY = {-Qm\ml - 20m|mt - Qm^ + m| s + 2ml m^s + ^^^t-^ + 4m|m^t 

/I 00 /I O O O 00 

-\-Amy^t + Sm^^m^u + 8my^,u — 2m^^su — 2m^su — Am^tu — 2m^u 
+su')/{2mUml-u)) 
TiuxTv = (64m|^m^m^ + QAm\m%ml - 16m|m^m|s - 32m|mtm^s 

— 16m^m|s + Am^jnlm^s^ + Am^m^s^ — 2m|^m^m|t — 28m|^m^m|t 

^10 /10 000 /10 /loo 

—2m^m2t + m^^m^st + 2m^^m^m2st + bm^m^st + 2m^^m^m2U 
—SGm^^m'^ym^^u + 2m^m|ii — mj^m^^su + lAml^mlm^^su + llm^m|sw 
—Amlmls'^u + Amf^mlm'^tu + Am'^m'^tu — 2mljn\stu — 2mlm\stu 
—Amj^mlm^u'^ — Am^m^u^ + 2m|^m^s-u^ + 2mlm\su^ — 2m^m|tM^ 

O*-? ^\// /I 0/0 \/0 \\ 

-\-m2:Stu + 2m^m2:U — m^su )/(2m^m2(m2 — s)[m{^^ — u)) 
TiiixTyi = (-64m|m^ - GAml^ml + 16m|^m^s + 32m|mts + 16m^s - Am\m1s^ 
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— 4m^s^ + 2mj^m'^t + 28m|^m^t + 2m^t — mj^st — — dm'^st 

—2m\^m^u + SGmf^m'lyU — 2m%u + mj^su — 14m|^m^sM — llm^sw 
+Aml^s^u — Amj^mlytu — Arnl^,tu + 2m'j^stu + 2m\,stu + Am^^jn^u^ 
-\-Am^u^ — 2rnj^su^ — 2rn^su^ + 2m^tu^ — stu^ — 2m^u^ + su^) 
/{2m%s{m\ - u)) 
= {l2m1-Amls + s'')/{Aml{m]j-s)) 
TJixTJv = (12mt-4m^s + s2)/(4mt(m^-s)) 

TivxTv = {—12m^m1.t + m^s^t + 12m^m|w — m^s^ii) / {Am^rn?^{m\ — s)) 
TiyxTyi = {12mit - sH - 12miu + s\i)/{4rnls) 

\T\^ = f^%x% + flTiixTii + flTnixTiii + flTiyX%j + /iTyxTy + fl%ixTyi 
+2fif2TixTu + 2/i/37ix7i„ + 2fjaixTiY + 2hU%x% + 2MsTuxTju 
+2/2/47JIX7JV + 2f2hriixTy + 2/3/47J„x7Jv + 2/3/57J„xrv 
+2UhT,yxTy + 2/i/67IxTvi + 2f2feTnxry, + 2MeTnixry, 
+2/4/6^vxrvi + 2MeTvxTvi (B2) 

Zl/J > 

Besides the £i exchange (III and IV in [17]), we also have 4 exchange, written here as VI 
and VII. 

VI. i-channel exchange 

VII. -u-channel £i exchange 

The couplings /i, . . . /s are modified, and we have new /e and /y couplings. 

/i = {-g2mzCos{/3 - a)/ cos9w){g2mz/ cos9w{{l/2 - sin^ 6w)cos{/3 + a) 

cos^ 9i + sin^ 9w cos(/5 + a) sin^ 9^) + g2m'j / {rriw cos /9) (— cos a) 

—g2'>Tie/ cos /?)(— ^4^ cos a — sin a) sin cos 
/2 = (-g'2"T.2sin(/5 - Q;)/cos^vt')(9'2"7.z/cos^ty((-l/2 + sin^ ^vk) sin(/3 + a) 

cos^ — sin^ 9^ sin(/9 + a) sin^ + (72"^^ / {^w cos /5) sin a 

—g2'>Tie/ {niw cos /?) (A^ sin a — /J, cos o;) sin cos 9^) 
h = (-Wcos^^(sin2^^-cos2^,/2))2 
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/4 = {-g2/coseyv{s\v?ew-cos'^9i/2)f 
/5 = -2gl/ cos^d^{sm^d^ + cos'^di{l/ A- sin^ 9^)) 
h = {-{g2/ cos9w)cos9tsm9i/2f 
h = {-{92/ cos9w)cos9iSm9i/2f 
TyixTyi = (m|^ - 4m|m| + Qmj^m'^z - Amj^ml + m| - 4m|^t + Am-mlt 

+Am\mlt - Amlt + 6m|t^ + Am\mlt^ + Gm^t^ - 4m|t^ - 4m|t3 

+i^)/(mt(m|-^)^) 
TviixTvii = (m|^ — Annq^m\ + &mqm\ — Arnj^m^^ + m\ — Am^^u + Awqw^^u 

+4m'j^m'^u — Am^^u + Gmj^u^ + Awq^m^v? + Qm\.v? — imj^u^ — Amlu^ 

W)/K(m| - 

TvixTvii = (m|^ + 12m|^m| + 38m|^m^ + 12m|^m^ + m| — 4:mj^m1s — 24m|^m^s 
-Amis + Amls"^ - '2m\t - lAm~m\t - lAm\m\t - 2m|i + Am\m\st 
+Am\st + m|^t^ + 2m|^m^t^ + m^t^ — 2m|^M — lAmj^mlu — lAmj^mlu 
—2m2U + Am^^m^su + Am^su + Am^tu + IGm^^m^tu + Am^tu 
—Am\stu — — + irij^u^ + 2m|^m|w^ + m^u^ — 2inj^tu^ 

-2mltu^ + iV)/(mt(m| - t)(m| - u)) 

. . / /10 0/1 (-1/1 00 /I 00 

TixTyi = [—6m^^m^ — 20mji^m^ — 6m^ + mji^s + 2m^^m^s + 5m^s + Sm^^m^t 

+8m'lt — 2mj^st — 2m|si — 2m\t^ + st^ + Arnj^rnlu + Am'^u — Am^tu) 

/(2mi{ml-s)(ml-t)) 
TiixTvi = (— 6m|^m^ — 20m|^m^ — 6m| + m|^s + 2m|^m^s + 5m^s + 8m|^m|t 

+8rr4t - 2m|st - 2m|st - 2m|t^ + st^ + Am]m\u + 4m^M - Am^tu) 

/{2mt{ml-s){ml-t)) 
TJnxTvi = -((m|^ + (m| - - 2m|(m| + i))7(mt(m| - t)(-m| + i))) 
TivxTvi = + + (6m| — t — m) + — m\{2s + t + u)Y 

l{m\{m\ -t){m\ - u)) 
TyxTvi = (— 6m|^m| — 20m|^m^ — 6m| + m|^s + 2m|^m|s + Sm^s + 8m|^m|i 

^%m\t — 2m| si — 2m|si — 2n?J^ + si^ + 4m| m^ii + Am\u — Am^^tu) 



41 



/(2mt(m|-t)) 



TixTvii = {—6m^^m^ — 20m^^mr, — Qmr, + m^^s + 2m^^m2S + bm^s + 4:m^^m^t 

A 00 A O O O 000\ 

+4:771 r^t + Sm^^m^u + 8m — 2m^^su — 2m ^su — Am^tu — 2m + su ) 
l{2m\{m\- s){m\-u)) 
TiixTvii = (— 6m|^m| — 20m|^m^ — 6m| + m|^s + 2m|^m|s + 5m^s + 4m|^m|i 

/I 00 A O O O 000\ 

+4m^t + Sm^^m^u + 8m ^u — 2m^^su — 2m ^su — Am^tu — 2m ^u + su ) 

/ {2m\{ml - s){m\ - u)) 
^iixTVii = (m|^ + + m| (6m^ - t- u)+tu- m\{2s + i + w))^ 

/(m^(m| -t){m\ - u)) 
TwxTyu = -((m| + {ml - uf - 2m|(m^ + u)f/{m\{m\^ - u){-m\^ + u))) 
TyxTvii = (— 6m|^m| — 20m|^m^ — 6m| + m|^s + 2m|^m|s + Sm^s + 4m|^m|i 

/I OO /t O O O OO 0\ 

+Am^t + Sm^^m^u + Sm^u — 2m^^su — 2m ^su — Am^tu — 2m ^u + su ) 
/{2m\{ml-u)) 

\rf - /i'7Jx7i + /|7iix7ii + /|7i„x7i„ + /27ivx7iv + /5Tvxrv + /6%ixrvi 
+/2rv„xrv„ + 2/i/27Jx7ii + 2/i/37ix7J„ + 2fjaixrw + 2f^f5TixTy 

+2/i/67Jxrvi + 2fJrTiXTyu + 2f2hTuy.%u + 2f2hTuXTiy + 2f2f5TiiXTy 

+2f2feTuxTyj + 2M7TuxTyu + 2hUTuixTiy + 2M,ruixry 
+2/3/6^iix7Vi + 2hhTnixTyn + 2UfaiYxTy + 2UfaivxTyi 
+2/4/7^vxrv„ + 2f,feTyxTyi + 2f,frTyxTyn (B3) 

There is no change for this channel, except that Tr — > £i. 
Zi^l — ^ ^7 

There is no change for this channel, except that Tr — > ^i. 
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ht^ — > Zh[H] 

Besides the exchange (I and II in [17]), we also have 4 exchanged, written here as IV and 
V. 

IV. i-channel 4 exchange 

V. li-channel 4 exchange 
The couplings are modified. 

fi = (-g2/cosey^{sm'^e^-cos^ei/2))(g2mz/cose^^((-l/2 + sm^ey^) 
sin[— cos] {P + a) cos^ 9e — sin^ sin[— cos] (/3 + a) sin^ 9e) 
+g2Tn'j/{mvi.' cos P) sin[— cos]q; — g2'mi/{rnw cos (3){Ae, sin[— cos]q; 
—/X cos[sin]Q;) sin Og, cos 6() 
f2 = -{-92/ cos6'vK(sin^ 9w - cos^ Oi/2)){g2mz/ cos9w{{-l/2 + sin^ 61,^) 
sin[— cos](/3 + a) cos^ — s\t? 6^ sin[— cos](/? + a) sin^ 
+g2'^1l iji^w cos /3) sin[— cos]q; — g2'me./{'m^ cos sin[— cos]a 
—^l cos[sin]Q;) sin 0^ cos 0() 
h = (-52/ cos0w{sm^ Ow - cos^6l£/2))(-g'2mzSin[cos](/3 - a)/ cosOw) 
fi = (-5'2Cos6'£sin6'£/(2cos6'vi.))(5'2^iz/cos^,,,((-l/2 + sin^6'v^,) 

sin[— cos] (/3 + a) (— cos Oi sin 0^) — siv? 9^ sin[— cos] {(3 + a) cos 9i sin 6'^) 
—g2'mi/{2mw cosP){Ae sin[— cos]q; — //cos[sin]Q;) cos(2^^)) 
fb = (-)(-5r2Cos6l^sin6l£/(2cos6lvK))(5'2"^z/cos6lvK((-l/2 + sin^6lvi,) 

sin[— cos](/5 + a){— cosOi sin ^^^) — sin^ 9w sin[— cos](/3 + a) cos9ism9i) 
~g2fnt/ (2mvK cos (3) {A^ sin[— cos]q; — /i cos[sin]Q;) cos(26'^)) 
TivxTJv = (ml + (ml-t)^ -2ml(ml + t))/{ml(ml-t)^) 
TyxTy = (m| + (m| - m)^ - 2m| (m| + M))/(m|(m| - m)^) 
TixTIv = -((m|^ + (m2 -t)2-2m|(m2+t))/(m2(m|-t)(-m|+t))) 
TiXly — (m|^ + + (6m| — i — + ixi — m|(2s + i + xi)) 
/(m^(m| -t)(m| - w)) 
TiixTiv = (m| + + m| (6m| — t — m) + — m^(2s + t + m)) 
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/(m^(m| -i)(m| - u)) 
TiixTv = -((m|^ + (m| - m)^ - 2m|(m| + M))/(m|(m| - M)(-m| + m))) 
TiuxTiY = {t{t - u) + mj^{-8ml - t + u) + 7nl{2s - t + u)) 

/{2ml{ml-s){ml-t)) 
TiuxTy = ((t - u)u + mj^{8ml - t + u) + ml{-2s - t + u)) 

/(2m|(m| -s)(m| - u)) 
TivxTv = (m|^ + + (6771! — i — m) + — m|(2s + i + «)) 

/(m|(m| -i)(m| - w)) 

+2fj2TixTu + 2fJsTixTni + 2hf^TixTiy + 2fJ,TixTy 
+2/2/3^ix7i„ + 2/2/47iix7iv + 2f2hTnxTy + 2f3famX%y 
+2/3/5^iixrv + 2UfaiYxry (B4) 

^ 7/l[if] 

There is no new channels, however the couphngs are modified. 

/i = (e)(5'2^z/ cos^^,4/((— 1/2 + sin^ ^^,4/) sin[— cos](/9 + a) cos^ 9e 

— sin^ 6*^. sin[— cos](/3 + a) sin^ 6'^) + g2'^'tl{fTiw cos /3) sin[— cosja 
—92'mil {rriw cos sin[— cosjo; — /i cosfsinjo;) sin 9^ cos 9e) 

f2 — — (e)(g'2"^2/cos^vK((— l/2 + sin^^vi')sin[— cos](/3 + q;)cos^^^ 

— sin^ 9w sin[— cos](/9 + a) sin^ 9e) + g2iTij/ {m^ cos /3) sin[— cosjo; 
—g'2?7i^/(?7ivv' cos /3)(74^ sin[— cosjo; — /xcos[sinjo) sin cos (B5) 

— > ZA 

We now have t- and li-channel ^2 exchanges, written here as III and IV. 

III. i-channel £2 exchange 

IV. M-channel £2 exchange 

The couplings /i and /2 are also modified. 

/i = {gifrizl cos 1/2 + sin^ ^^y) sin(/3 + o) cos^ — sin^ 9yv sin(/3 + o) sin^ 9i) 
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+g2m'j I {m^^ cos /3) sin a — g2Tni/ {m-w cos (5){A(^ sin a — /i cos a) sin $£ cos 
(—92 cos(a — /5)/(2cos^w')) 
/2 — (5'2^z/ cos6'vv.((l/2 — sin^ 9^) cos(/3 + a) cos^ 9i + sin^ cos(/3 + a) sin^ 6*^) 
—g2'>Ti£ I {vrtyf cos /3) cos a — g2fn(,/ {ruw cos (3) {—A^ cos a — sin a) sin 9^ cos 
(-g'2 sin(Q; - (3)/ (2 cos ^vk)) 
/s = (-5'2Cos6'^sin6'|./(2cos6'vi.))(-5'2"^£/(2miv)(A^ tan/3 - /x)) 
/4 = {-g2Cos9ism9£/{2cos9yv)){-g2rn£/{2rn^){A£iaj\(5 - jj)) 
TuixTui = (m| + (ml - tf - 2m|(m^ + i))/(m^(m| - i)') 
TiyxTiY = {mj^ + {ml - uf - 2m|(m| + (m|(m| - uf) 
TmX%y = (m|^ + m'l + (6m^ — t — m) + — m|(2s + t + «)) 
/(m|(m| -i)(m| -u)) 
TjxTin = (— + + m|^(— 3m^ + s) + — si + m^(— — 3m| + 

+2m|M)/(m|(m^ - s)(m| - t)) 
%iX%\i = (— + mis + (— 3m| + s) + m^t — st + m^(— — 3m| + t) 

+2mlu)/{ml{m'jj - s)(m| - t)) 
TJxTiv = (~''^z + "^z'^ + mj^{—3ml + s) + 2m|i + — + m\{—m'j^ — 3ml 

+u))/{ml{ml - s)(m| -m)) 
TiixTiY — {—m'l + mls + m'j^{—3ml + s) + 2mlt + mlu — su + ml{—mj^—3ml 
+u))/{ml{mjj - s)(m| - u)) 
|r|2 = f^T.xT, + fiTuxTu + /iTinxTIn + f!T,yxT,y + 2fj2TjxTu + 2fJ,TjxTui 

+2f,UTixTiy + 2f2hTiixTiii + 2/2/4T„x7iv + 2hf^TuiX%, (B6) 

^1 > 

The s-channel /i and if annihilations are neglected, due to the small Yukawa couplings for 
leptons. The ^-channel x exchange is modified, with more couplings introduced {K, K') — > 
[Ka, Ki,, K'^, Kl). The terms with m^ were neglected in [17], but not here. The couphngs fsc 
and fad are modified while f^c remains the same. 

/sc = {-g2{sm'^9^ - cos^ 9^/2)/ cos9^){g2{l - ^siv? 9y^)/{^cos9^)) 
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fsd = (-^2(sin2 0^ _ ^Qg2 cos0w){-g2/iAcose^)) 
Ka = -s\\i6ii{g2rni>/(2mwCosf5)Ni:>, + giNii)/\/2 

- cos 6e{g2mi/ {2m^ cos P)Nis - {giNn + g2/2{Ni2 - tan^^)))/x/2 
Kb ^ - sin Oe{g2mi/ {2mw cos /3) A^is - 9iNii)/V2 

- cos Oi{—g2mi/ {2mw cos 

= -sin^^(^2"^£/(2m^cos/3)iVj-3 + ^iiVii)/V2 

- cos Oi{g2mi/ {2mv^ cos 

= -sin^^(^2m^/(2mv^cos/3)^,-3-^i%)/v^ 

-cos^f(-^2m^/(2mv^cos/3)iV,-3 - (^iiV,i + ^2/2(iV,-2 - iV,i tan ^>^)))/a/2 
TJiixTV = {2f^dKbKa{-lQm]m] + 4m|s + 4m^s - + - 2tM + u^) 

+/3c(-f^a (~4m^^s + s — Arrig^t — Amem^qt — t + Am^u + Amem^qu + 2tu — u ) 
-\-Kl{—Amj^s + — Amit + Am^m^qt — t^ + Amju — Am(,m^qu + 2tu — u^))) 
/{{ml-s){m\,-t)) 

Tj^xTy — —{f4c{Kl{—Amf^s + s'^ — Am^t — Amim^qt — t'^ + Arn'ju + Arni'm^qu + 2tu 
—u ) + Kf^ {—Am^^s + s — Am^t + Amem^qt — t + Am^u — Amim^qu + 2tu 
-u'))/{s{mlo-t))) 

TyXly = {2{K'^K^{—mj^ — 2m|^m^ + 3m^ — mfm^o — Smfm^q + mfm^q + 3mfm^q 
—2rrijm^qm^q — 2rrijm^qm^q + rnj^s — mjs + mim^qs — mim^qs 
-\-m^qm^qs + m-^t + 2m^t + m^t — mem^qt + mem^qt + m^^u — m^u 
+mim^ou — mim^ou — tu) + K^KaK^Kb{—2'm^^ — Am^^m^ + 6m^ + 2m^^s 
—2m(S + 2m^qm^os + 2m^^t — Am^t + 2'm^t + 2m^^u — 2m^u — 2tu) 
{—m^^ — 2m^^m£ + Sm^ + Am^m^q — Am^m^q — Am^m^qm^q 
+m|^s — mjs — mj^m^qs + m^m^qs + m^om^os + mjt + 3m|i + m(^m^qt 
—m^m^ot + m|^M — m^w — rriim^ou + mem^ou — tu) 
+ {KaK + KbKl){KaK'^ + KbK'^){-m]s - m^om^os + mj t + mjt 
+mju - mju - tu) + {KaK + KbK'^){KaK - KbK'^){Am\m^.^ 
+Am^m^q — mgm^os — mim^os + m^m^ot + m^m^qt — me,m^qu — mim^qu) 
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/((m|o-t)(-m|o+t)) 
|r|2 = 7i„x7i„ + 7ivx7iv + Ervxrv + 27i„x7iv + 2^(7i„xrv + 7ivxrv) (B7) 

Ul ff 

When / is a fermion other than a lepton or t quark, again the s-channel h and H annihilations 
can be neglected. We separate the channel Iill — > viVf^ because it has i-channel chargino 
exchange. The couplings jzc and /s^ are modified while /4c remains the same. 

he = (-y2/cos^^(sin2^^-cos2^,/2))(^2(-2T3^ + 4g_^sin2^^)/(4cos^^)) 

/3d = (-^2/cos^^(sin2^^-cos2^,/2))(y2(2T/)/(4cos^^)) (B8) 

lxl\ > ViVi 

With i-n wc have t- channel chargino exchange. 

III. s-channel Z annihilation 

IV. t-channel X{\ 2) exchange 

/s = (-5'2/cos6'^(sin26'i^ - cos^6'^/2))(-5(2/(2cos6'm.)) 
/4(i) = {g2mi/ {V2mw cos (3)Ui2sm 9f ~ g2Uii cos 9 
TiuxTui = (-4m|s + + 2tu -t^ - v?) /{s - mlf 

TnixTiy = -{-2m}s + l/2s^ + tu-l/2t^-l/2u^)/{{s-ml){t-ml+)) 

1 /Ci 

\T\^ = fiTuixT,u + Y.f4ii)U{j)TjvxTiy + 2hY,Uii)TuixTiy (B9) 



ilil ^ tt 

The couplings fia, /2a, fsc and /3d are modified while /4c remains the same. 

/la = (— 5'2"^t sin a/(2miv sin/9))(5f2^z((l/2 — sin^ cos(/3 + a) cos^ 
+ sin^ 9w cos{P + a) sin^ 9£) / cos 9w — 5'2?77.| cos Q;/(mvv. cos (5) 
—g2mt{—A^ cos a — sin a) sin cos 9^/ {rriy^^ cos pi)) 
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f2a = {-g2mtCOsa/{2'mwSmP)){g2rnz{{-l/2 + sm'^9w)sin{P + a)cos^9e 
— sin^ 9w sin(/5 + a) sin^ 9i) j cos 9w + Qif^t ^wiajimw cos /?) 
~5'2^£(^£ sin a — /J, cos a) sin 9^ cos 6*^/ (mv^ cos 
/sc = (-^2(sin'^^-cos'^^/2)/cosev^)(^2(-l + 8/3sin'e^)/(4cosev^)) 
/3d = (-y2(sin'^^-cos2^,/2)/cosM(W(4cosM) (BIO) 

iJl > hh 

The couplings are modified. Channels IV (t-channel £ exchange) and V (w-channel £ ex- 
change) are now summed over £i and £2, with appropriate propagators. 

/i = (-35f2"^z cos(2q;) sin(Q; + /9)/(2cos6'vk))(5'2'tT'z((1/2 - sin^^iv) cos(/5 + a) 
cos^ 9e + sin^ 9^ cos{P + a) sin^ 9i) / cos 9-^^ — g2'iTij cosa/{m^ cos (5) 
—g2fni{—'^t cos a — II sin a) sin cos 9^/ {■my^r cos /5)) 
/2 = ((72'^z(cos(2a) cos(a + /9) — 2 sin(2a) sin(a + /5))/(2 cos^^iv)) 

(5'2^z((l/2 — sin^ 9w) cos(/3 + a) cos^ 9e + sin^ cos(/3 + a) sin^ 6*^)/ cos6'vk 
—g2Tnl cos Q;/(mvK cos /3) — g2mt{—A(^ cos a — // sin a) sin 9^ cos 9(,/{myv cos /3)) 
/3 = (t/^/2)((-l/2 + sin^^H/)cos(2Q;)/cos^^v^ -m^sin^Q;/(m^cos^/9))cos^^^ 
-|-((72/2)(— sin^ 9w cos(2q;)/ cos^ — m\ sin^ Q;/(m^ cos^/5)) sin^ 9^ 
/4(1) = (gi2m2((l/2 — sin^ 9yv) cos{P + a) cos^ 9(, + sin^ cos(/9 + a) sin^ 6'^)/ cos6'vk 

—g2'>Tij cosa/{my^^ cos /3) — g2m(,{—A(^ cos a — // sin a) sin 9^ cos 9(,/{my^^ cos /3))^ 
/4(2) = (5(2m^cos^£sin^^((l/2 — sin^^^^) sin(/9 + a) — sin^ sin(/3 + a))/ cos9w 

—g2^i{Ai sin a — yucoso;) cos(2^^^)/(2m,^ cos/3))^ 
/5(1) = {g2mz{{l/2 — sin^ ^vv) cos(/9 + a) cos^ + sin^ 9w cos(/? + a) sin^ 9e)/ cos^vv^ 

—g2Tnl cos Q;/(mvK cos (5) — g2m(,{—A^ cos a — // sin a) sin cos 9£/{myv cos /3))^ 
/5(2) = {g2mz cos ^£sin^£((l/2 — sin^ 9w) sin(/9 + a) — sin^ ^^^v sin(/3 + a))/ cos6'v^ 
—512771^(74^ sin a — n cos a) cos(26'£) / {2m^v cos /3))^ 
|r|^ = f^TixTj + f^TuxTu + fiTuixTui + Y.{Mi)f4{j)riyxTjy + Mi)M^^^ 

+2fj2TixTu + 2fJ,T,xTni + 2^(/i/4(7)7ix7iv + hh{^r^y-%f) 

i 
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i 

+2 E(/3/4(^)^iiX^v + hh{^ruix%j) + 2 5: h{i)h{3)%,x% (Bll) 
— > hA[HA] 

Besides the s-channel Z annihilation, we have t- and w-channel I2 exchanges. 

II. i-channel £2 exchange 

III. li-channel £2 exchange 

The couphng /i is also modified. 

/i = (—5(2 cos[sin](a — /5)/(2 cos^^i4/))(—5'2(sin^ 6^14/ — cos^ ^^£/2)/ cos^^w.) 

h = (-5'2m^(^£tan/?-//)/(2m;v))(c/2mz((l/2-sin^6'iy)sin[-cos](/? + Q;) 
— sin^ 9^ sin[— cos](/3 + a)) cos sin cos 9^ 
—g2'>Tie{Ae sin[— cos]q; — fi cos[sin]Q;) cos{29e)/ {2mn' cos jS)) 

/s = (-5'2?^i£(A/-tan/5 - /i)/(2m^^,))(5(2m2((l/2 - sin^6'vy) sin[-cos](/? + a) 





— sin^ 9w sin[— cos](/3 + a)) cos sin cos^vv 




—g2me{A£ sin[— cos]q; — cosfsinja) cos(2^^) / {2mw cos 




= i/K-^r 


TuixTui 


= l/(m|-«r 


TixTii 


= (i-ti)/((m^-s)(m|-i)) 




= (i-'u)/((m|-s)(m| -«)) 


TuxTui 


= l/((m|-t)(m|-«)) 




= /^TJxTJ + /|r„x7ii + fiTuixTui + 2hf2%x%^ + 2hU%xrni 




+2f2f3TuxTui 



ilil — ^ W+H- 

Besides the s-channel H and h annihilation, we have t-channel sneutrino exchange. 

111. t- channel z/^ exchange. 

The couplings /i and /2 are modified. 

fi — {g2iTiz{{l/2 — sin^ 9w) cos(/5 + a) cos^ 9i + sin^ 9w cos{P + a) sin^ 9i)/ cos9. 
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—g2m1 cosa/{mw cos /3) — g2mi{—A(, cos o; — sin a) sin 6(, cos (mvv^ cos P)) 
(^2sin(«-/3)/2) 

/2 = (S'2^z((— 1/2 + sin^ 9y^) sin(/3 + a) cos^ 6*^ — sin^ 9^ sin(/3 + a) sin^ 6()/ cosO^ 
+g2'>Ti'j sm a /{rriw cos /3) — g2mt{At sin a — ji cos a) sin cos (mvv^ cos 
(^2Cos(a-/3)/2) 

h = (-5'2/\/2cos6'<?)(-5f2"^vi'Cos6'£(-m^tan/3/m^ + sin(2/3))/\/2 
+92me/ {■\/2m^y) sin 6'^(/^ - tan/3)) 
ruixTui = (m|^ + (m^ - tf - 2m\ {ml + t))/{ml{m% - tf) 
TixTiii = {-ml + m^s + m| (-3m^ + s) + m^t - si + m|+ - 3m^ + t) 

+2mlu)/{ml,{mjj - s){ml - t)) 
^ixTJii = {-ml + mis + m| {-3ml + s) + mlt - st + m]j+ (-m| - 3m^ + t) 
+2mlu)/{ml{ml - s)(m?^ - t)) 
ITf = f^TixTi + /2'7Iix7Ii + fiTuixTui + 2fj2TixTn + 2hhTixTni 

+2/2/3^ix7i„ (B13) 

The new channels are 

IV. t-channel I2 exchange 

V. M-channel £2 exchange 
The couphngs are modified. 

/i = {g2mz{{l/2 — sin^ 9^,) cos(/3 + a) cos^ 9e + sin^ cos(/3 + a) sin^ ^^) / cos 9w 

—g2mj cos a/ (m^y cos 0) — g2m^{—A£ cos a — sin a) sin6'f cos 9i/ {m^y cos /?)) 

(5'2"^zCOs(2/?) cos(q; + (3)/ {2 cos 9w)) 
f2 — {g2mz{{—l/2 + sin^ 9-^^) sin(/3 + a) cos^ — sin^ 9^ sin(/3 + a) sin^ 9i) / cos^vv 

-\-g2ml sin Q;/(miv cos 13) — g2mi{Ae sin a — (i cos a) sin^^^ cos 9i/ {m^ cos /3)) 

{-g2mz cos(2/3) sin(Q; + (3)/{2 cos ^i^)) 
/3 = ^^/2((-l/2 + sin^ 9w) cos(2/3)/ cos^ 9^ - m] tan^ (3 /ml) cos^ 9e 

+g2/2{- sin^ 9^ cos(2/3) / cos^ 6*,^ - m| tan^ /3/m^) sin^ 6*^ 
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/4 = {-g2me{Aiteinl3 - ii)/{2mnr)y 

fb = {- g2me{Ai tan P - ij,)/{2mv^))'^ 

T^vxTiy = l/(m|-t)2 

TvxTv = l/{ml-uf 

r.xTiy = i/((rnl-s)iml-t)) 

Tux%y = i/((rnl-s){ml-t)) 

TuxTy = i/{{rnl-s){ml-u)) 

TiuxTiy = l/(m|-t) 

TmxTv = l/(m|-M) 

|r|2 = f^T,xTj + f^TuXTu + fiTuiXTui + f!T,yXTjy + flTyXTy 

+2hf2%x%i + 2hfaix%ii + 2fjaixTjy + 2fjaixTy 

+2f2hTuxTjn + 2f2hTiixTiy + 2f2hT^ixTy + 2hUTuixT,y 

+2hhrmx%j + 2Uf,TiYxTy (B14) 

The and w-channels £ exchanges are now summed over Ii and with appropriate propa- 
gators. The couphngs are 

fi — {g2'iTiz{{^/2 — sin^ 9w) cos{P + a) cos^ 9i + sin^ 9^ cos(/9 + a) sin^ cos^^y 
—g2mi cosa/{m^ cos /3) — g2mi{—Ai cos o; — // sin a) sin 9^ cos (m^^. cos (3)) 
((?2^z(cos(2a) sin(a + /9) + 2 sin(2a) cos(« + l3))/{2 cosOw)) 

/2 = (5'2'^z((— 1/2 + sin^ 6*1^) sin(/5 + a) cos^ 6*^ — sin^ 9w sin(/3 + a) sin^ ^^^)/ cos6'w- 
+g2ml sin a/ (mp^. cos /3) — g2mj({Aj^ sin a — /i cos a) sin cos (m^^. cos (3)) 
{g2'mz{cos{2a) cos{a + P) —2 sin(2Q;) sin(Q; + /3))/(2 cos 9w)) 

fs = ^2sin(2a)((-l/2 + sin2^^)/cos2^^^ +m2/(2m^cos2/3))cos2^f/2 
+gl sin(2Q;)(- sin^ 9^/ cos^ 9^ + m^/(2m^ cos^ (3)) sin^ 61^/2 
/4(1) = (5i2?7iz((l/2 — sin^ 9w) cos(/3 + a) cos^ 9^ + sin^ ^p^. cos(/3 + a) sin^ cos^^y 
—g2m\ cosa/{myv cos /3) — g2m£{—A(^ cos o; — // sin a) sin 9(^ cos (m^^. cos (3)) 
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(fl'2"^z((— 1/2 + sin^ 9w) sm(/3 + a) cos^ 9e — sin^ 9w sin(/3 + a) sin^ 9()/ cos9w 
+g2mj sin a/ (ruw cos /3) — 5f2^£(^^ sina — fi cos a) sin cos 9£/ {rriw cos 
/4(2) = (g(2mz((— 1/2 + sin^ 9^) cos(/3 + a) + sin^ cos{j3 + a)) cos 6*^ sin 9^/ cos9^y 
—g2mi{—Ae cos a — sin a) cos(2^^) / (2mvi. cos 

(fl'2"T.z((l/2 — sin^ sin(/9 + a) — sin^ sin(/9 + a)) cos sin ^^/ cos^v^. 
— fl'2'^^(^^ sina — yucoso;) cos(2^^£)/(2m^y cos/3)) 

/5(1) = {g2mz{{l/2 — sin^ 9yy) cos{P + a) cos^ 9^ + sin^ cos(/3 + a) sin^ ^^)/ cos^vv' 
—g2Tn1 cosa/im^ cos /3) — g2m(^{—A^ cos o; — // sin a) sin 9^ cos ^^/ (mvv. cos (3)) 
(fi'2"^z((— 1/2 + sin^ ^vi') sin(/9 + a) cos^ — sin^ 9w sin(/9 + a) sin^ ^f) / cos^v^ 
+g2m'j sma/{myy cos /?) — g2m£{A£ sin a — /i cos a) sin cos9e/{my^, cos /?)) 

/5(2) = {g2mz{{—l/2 + sin^ ^vk) cos(/9 + a) + sin^ cos(/9 + a)) cos9esm9e/ cos 9w 
—g2me{—A£ cos a — //sin a) cos(2^^)/(2m^y cos/3)) 

(5'2"^z((l/2 — sin^ 9w) sin(/? + a) — sin^ 9w sin(/3 + a)) cos 9e sin ^£/ cos^vv 
—g2mi{A£sma — /^cosq;) cos(26'£)/(2mvK cos/3)) 
\T\' = f^TjxTj + f^TuxTu + fiTuixTui + Y.{Mi)f4{j)riyxTjy + Mi)M 
+2fj2TixTn + 2fJ,T,xTni + 2J2 f,{U{t)T,xT,y + Mi)T,xTy) 

i 

+2f2hTnX%u + 2 5^ f2{f4{i)ruxTjy + Mi)TnxTy) 

i 

+2/3E(/4(0^IlX^V + hmuiXTy) + 2^U(l)Mj)riyXTy (B15) 

The t- and w-channels £ exchanges are now summed over ii and £2 with appropriate propa- 
gators. The couphngs are 

fi — {—3g2mz cos{2a) cos(q; + /3)/(2cos^vv'))(fl'2"^z((l/2 — sin^ 9w) cos(/3 + a) cos^ 9e 
+ sin^ 9w cos{P + a) sin^ 9i) / cos 6*^^ — g2'm'j cos a /{ruw cos /3) 
—g2'm'e{~^i cos a — II sin a) sin cos ^^/ (mv^. cos /?)) 

/2 = (g'2"^z(cos(2Q;) sin(Q; + /3) + 2sin(2Q;) cos(q; + /3))/(2cos^vk)) 

{g2'mz{{—l/2 + sin^ ^vv) sin(/? + a) cos^ 9e — sin^ 9wSm{P + a) sin^ 9e)/ cos9w 



+g2m'j sma/{'mw cos (3) — g2mi{A£ sin a — /i cos a) sin 9^ cos 9i/{mw cos (5)) 
/a = (?2/2(—(— 1/2 + sin^ ^iv) cos(2a)/ cos^ ^v^. — cos^ o;/(m^ cos^ cos^ 

+gi2/2(sin^ 6*^^ cos(2q;)/(cos^ Oy^) — 'mj{cos^ a)/{rn^ cos^ sin^ 6*^ 
/4(1) = (gi2"T.z((l/2 — sin^ ^vk) cos(/3 + a) cos^ + sin^ 9^ cos(/3 + a) s\v? 9i)/ cos9w 

—g2m'j cosa/{mw cos /3) — g'2"^£(— cos a — // sin a) sin cos 9i/{m^ cos 
/4(2) = (5(2777.2 ((1/2 — sin^ cos(/3 + a)(— cos sin 6^) + sin^ ^^t^ cos(/9 + a) cos 9i sin 6'£) 

/ cos 9^ — g2mi{—A( cos a — ji sin a) cos(26'^) / (277ivv' cos (5)Y 
75(1) = {g2mz{{l/2 — sin^ ^vk) cos(/3 + a) cos^ + sin^ 9wC0s{P + a) sin^ 9i)/ cos 9yy 

—g2m'j cos a /{myy cos (3) — g2mi{—Ai cos a — fisin a) sin 9i cos 9i/{mw cos P)Y 
75(2) = {g2mz{{l/2 — sin^ 9^) cos{(3 + «)(— cos 6*^ sin 6'£) + sin^ 9^, cos(/5 + a) cos 6*^ sin 6'^) 

/ cos 9yy — g2mi{—Ai cos a — // sin a) cos(2^^) / (2777vv. cos /?))^ 
|r|^ = f^%xTj + f^TuxTu + fiTuixTui + J2(Mi)Mj)'rwxTjy + Mi)M^^^ 

+2fif2TixTu + 2/i/37ix7i„ + 2^(/i/4(z)7ix7iv + fiM^T.xTy) 

i 

+2f2hTuxTni + 2 Y.ifM^'^nX'^iy + f2h{i)riixTv) 

i 

+2Y,{hfA{i)%iiX%j + hh{i)%iix%) + 2Y,f4{i)Mj)riyxTy (B16) 

— ^ H+H- 

We have one new channel, 

VI. f-channel exchange 

The couplings, except f^, are modified 

/i — {g2fTiz{{l/2 — sin^ 9w) cos(/3 + a) cos^ 9i + sin^ 9w cos(/3 + a) sin^ 9()/ cos9w 
—g2'>Ti'j cosa/{m.yi, cos (5) — g2m£{—A£ cos a — /i sin a) sin 9^ cos 9(,/ {m^ cos (5)) 
{—g2{mw cos{P — a) — ruz cos{2P) cos(/3 + q;)/(2cos^vk))) 

/2 = (5'2'^z((— 1/2 + sin^ 9w) sin(P + a) cos^ 9£ — sin^ 9w sin(/9 + a) sin^ 9i)/ cos9w 
+g2m£ sin a/{mw cos P) — g2mi{AiSina — //cos a) sin cos ^^/ (777 vk cos/?)) 
{-g2{mw sin(/3 - a) + 777^/ (2 cos 9yy) cos(2/3) sin{(5 + a))) 

/a = (-g'2(sin^6'iv - cos^6'£/2)/cos6'tt.)(5(2Cos(26'vt^)/(2cos6';^)) 
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/s = gl cos(2/3)(l + (-1 + 2 sin^ 9^)/ {2 cos^ 9^)) cos^ 9^/2 

-{gl cos(2/5) sin^ 9^/2 + glm] tan^ /?/(2m^)) sin^ 
/e = (-5f2"^w' cos 6'^(-m^ tan /3/m^ + sin(2/?))/\/2 
+g2mesin9e{n - A^teni f3)/{V2mw)y 
TvixTvi = l/(m?^-t)' 
TIxTvi = l/((m^-s)(m2^-t)) 
TiixTvi = l/{{ml-t){ml-t)) 
^iixTvi = (t-ii)/((m^-s)(m?^-i)) 
TJvxTvi = {-t + u)/{s{ml-t)) 
TvxTvi = l/(m?^-t) 

+2/1/27^x711 + 2/i/37Jx7J„ + 2/i/47Jx7Jv + 2hfaiX% + 2hfaiX%i 

+V2hTuXTui + 2hhTllX%y + 2hhTuxTy + 2hhTuXTyi 

+2f3hTmxTiv + 2/3/57i„xrv + 2Mamy~%i + 2UfaiyX% 

+2/4/6^vXrvi + 2f,feTyXTyj (B17) 

^ 

The channels I. i-channel x exchange II. w-channel x exchange are modified as follows. 

= — sin^£(g'2"T.£/(2mM'COS 

— COS {g2'me / {2m^ cos /3)iV,3 + (-^liVa - ^2/2(iVi2 - iVatan^„,)))/y2 
Xfe = -sin6'^(c/2?nV(2mw'Cos/?)A^i3 - c/iA^a)/A/2 

— cos 9(,{—g2m(,/ {2m^ cos /3)A^«3 + {-giNi^ - g2/2{Ni2 - Nnte.n9^)))/V2 
< = -sin^^(^2"^^/(2m^cos/5)Ar^-3 + ^iiVii)/\/2 

— cos9e{g2me/{2m^ cosP)Njs + {-g^N^i - g2/2{Nj2 - Nj^ tan9^)))/V2 
K\j — — sin 9 e{g2me/ {2171^008 

— cos 9(,{—g2m(,/ {2m^ cos /3) AT.g + (-5,7V., - g2/2{N^2 - Nj, tan 9^)))/V2 
TixTi = {-2{AKaKKbKm^om^o,s 
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+m|^s — mjs + rriim^os — m^m^os + m^om^os + mqj^ + Sm^i — mim^ot 
+m£m^ot + m'j^u — m'ju + mem^ou — m^m^ou — tu) 

—nrj^s + mjs — rriim^qs — rriim^qs + m^qm^qs — mj^t — Smjt + rriim^qt 
+mem^ot — mj^u + mju — rriim^ou — m^m^ou + tu) 

{—'m^^ — 2m^^m^ + 3m^ + 4m^m^o — Am^m^q — Am^m^qm^q 
+m|^s — mfs — mfTn^qs + rrifrn^qs + m^om^gs + m|^i + Sm^i + m(,m^qt 
—m^m^qt + — — m^m^qu + m^m^qu — tu) 

(m^^ + 2m^^m^ — 3m^ — 4m^m^o — 4m^m^q — Am^m^om^q 
—inj^s + mis + rrif^m^qs + rrif^m^qs + m^^om^^gs — m|^i — Sm^i — m(^m^qt 
—m(,m^qt — mf^u + mju + m(,m^qu + m(,m^qu + iii))) 
/((m|o-t)(-m|o+t)) 

— 2mj^mj + 3mj + Am^m^q — Am^m^q — Arn^m^qm^q 
+m|^s — mjs — rriim^qs + rriim^qs + m^om^qs + mf^t — mjt — m(,m^qt 
+mim^ot + mj^u + Sm'ju + rriim^QU — mim^nu — tu) 

(m^^ + 2m^^m^ — 3m^ — 4m^m^o — Am^m^q — Am^m^qm^q 
— m|^s + + rrif^m^qs + rrif^m^qs + m^^gm^^gs — 77T.|^i + + rrif^m^qt 
+mimj^ot — mf^u — Smju — rriim^ou — rriim^qu + tu) 
+K^K'^{—mj^ — 2rrij^rnj + Smj — Am^m^q + Am\m^a — Amjm^qm^q 
+m|^s — mis + m^m^qs — m^m^qs + m^qm^qs + m'j^t — m^t + m^m^qt 
—m^m^qt + m^^u + 2>mlu — mim^qu + mim^qu — tu) 
+Kf^K,j [m^^ + 2m^^mi — 3m^ + Am^m^o + Am^m^o — Am^m^om^o 
—mj^s + mjs — m^m^qs — m^m^qs + m^qm^qs — rnjt + m^t — m^m^qt 
—m^m^qt — mf^u — Sm^u + me,m^qu + m^m^qu + tu))) 
/((m|o - ii)(-m|o + u)) 
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TjxTu = {SKaK'^Ki,K'f,m^mi^os + 

K^{—Kfj {2m^^ + Arriijn^ — 6m^ — d>m^m^o + Smim^o + d>'m^m^om^o 
+2m^^s + 2m^s + 2m^m^os — 2m^m^qs — 2m^om^qs — s — 2m^^t — 2mgt 
+2m^m^oi + 2mem^qt + t— 2m^^u — 2m^u — 2mem^qu — 2mem^qu + u )) + 
Kg^K^ {2m^^ + Am^m^ — 6m^ — Sm^m^q — Sm^m^q — %m^m^qm^q 
+2m^^s + 2m^s + 2m^m^os + 2m£m^os + 2m^om^qs — s — 2m^^t — 2m^t 
+2mim^qt — 2mim^qt -\-t^ — 2rrij^u — 2mlu — 2m(rn^qu + 2mi'm^qu + u^) + 
K^K^ {2m + Arrij^m^ — 6m ^ + Sm^m^q + Sm^m^q — Sm^m^qm^q 
+2m|^s + 2m'js — 2m(^m^qs — 2m^m^os + 2m^am^as — — 2m|^t — 2m|i 
—2mem^at + 2mem^ot + — 2m^ -u — 2m'ju + 2mfm^au — 2mg,m^ou + m^) — 

Ai tl ^ A.i A.j ^ 

K^K'^{2mj^ + iwqjnl — 6mj + 8m^m^q — Sm^m^q + Sm^m^qm^q 

o o o o o 

+2m^^s + 2m^s — 2mem^qs + 2mim^qs — 2m^om^qs — s — 2m^^t — 2m^i 

O O O O \ \ 

—2mim^ot — 2mim^)t + t— 2mg^u — 2m^u + 2mim^qu + 2mem^qu + u )) 
/((m|o-t)(m|o-«)) 

\T\' = J^(T,xT, + T,,xT,, + 2T,xTu) (B18) 

If If* > i^fB 

I. t-channel x exchange 

Here mA,B = m^A.s and m^^ = m^A,B. 

A{i) = -g2/V2{cos9A{-Ni2 -tan 9w Nil) + sin 9 Am a/ {mw cos P)Nis) 

B{i) = - g2 /\/2{cos 9 Am A /{m COS pi) A^j3 + 2 sin 9a tan 9^Nii ) 

C{i) = — g2/V2{cos 9 BmB/{mw cos (3) + 2 sin 9b tsm 6 w Nil) 

D{i) = - g2 / \/2{cos 9 b{- Ni2 - i&n 9^N ii) + sin 9 BmB/{m^ cos j3)Ni^) 
TixTi — {A{i)C{i)A{j)C{j){—m\myqm^q—m%myqm^q+myqm^qs) 
+A{i)C{i)A{j)D{j){—m'~mBm^q +m\mBm^o + mBm^ot) 
+A{i)C (i) B (j)C (j) {m~mAm^q + 2mAm^^m^q — mAm^qs — mAm^qu + m^m^o) 
+A{i)C{i)B{j)D{j){-2mAmBm^om^o) 
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+A{i)D{i)A{j)C{j){—'m'^^mBm^o -\-m\mBm^o -\-mBm^ot) 

A/\/\i/\/\/ OO O O O A O 

+A(i)D{i)A[j)D[j)[2rn~rnB — m~s — m~t — + m~ — 2mj<jnB + mj<^s 

+m\u — m\ — m\t + tu) 

+A{i)D{i)B{j)C{j){-2mAmBt) 

+A{i)D{i)B{j)D{j){rn~mAm^o + 2m^m^m^p — niAin^qs — niAm^qu + ni\m^q) 
+B{i)C{i)A{j)C{j){m'~mAm^o + 2mAm^B''^x° ~ TnAm^os — mAm^ou + m\m^o) 
+B{i)C{i)A{j)D{j){-2mAmBt) 

+B(i)C{i)B(j)C{j){2nt~m% — m~s — rrt~t — m-u + m~ — 2m^m^ + m\s 
-\-m\u — m\ — m\t + tu) 

+B{i)C{i)B{j)D{j){~m~rnBm^q + m\mBrn^o + niBm^qt) 
+B{i)D{i)A{j)C{j){—2mAmBm<^qm<^q) 

+S(i)£)(i)74(j')£)(j')(m~m^m^o + 2mATn\m^q — mAm^qs — mAm^qu + m^m^o) 
+B{i)D(i)B{j)C{j){—m'~mBm^o + m\mBm^o + mBm^ot) 
+B{i)D[i)B[j)D[j)[—m\m^qm^q — m\m^qm^q + m^qm^qs)) 
l{{t-m\o){t-m\o)) 

|r|2 = J^TJxTJ (B19) 

pApB pApB 

1 1 ^ 

I. t-channel x exchange 

A{i) = -c/2/ a/2(cos 9A{-Ni2 - tan Oy^Nn) + sin OAmA/ {m^ cos /3) A^js) 
i?(z) = — (72/v^(cos^yimyi/(miy cos/9)iVj3 + 2sin^^tan^^yA^ji) 

C(i) = — (72/A/2(cos^^BmB/(mvy cos/3)iVi3 + 2 sin^B tan^iv^ii) 
-D(i) = -g'2/A/2(cos6'B(-A^j2 -tan ^vv'-^ji) + sin 6'BmB/(mvK cos/?) A^is) 
TixTi = {-A{i){-B{j)mAiD{i)m^o{2C{j)mBm^o + D{j){m\ + 2m| + m~ - s - w)) 
+C(i)(2L)(j)mBt + C{j)m^o{m\ + 2m| + m| - s - «))) 
+A{j){D{i)rn^o{D{j)rn^o{rn\ + m% - s) + C{j)mB{m\ - m^~ + t)) 
+C {i){D{j)mBm^q{'m\ — m- -\-t) -\- C{j){m\ — m~ + m-s + 
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+m|(-2m~+ + m\{2m% - s - + m~ii - tu)))) 
-B{i){C{i)m^o{B{j){C{j)m^o{m\ + ml - s) + D{j)mB{m\ -m\ + t)) 
-A{j)mA{2D{j)mBm^o + C{j){m\ + 2m| + m^~- s- u))) 
+D(i)(C(j)mB(-2A{j)mAt + B{j)m^o{m\ - m\ + t)) 

+D{j){A{j)mAm^o{—m\ — 2m| - m~ + s + u) + B{j){m\ — m~ + m~s + m'-t 

+m|(-2m~ + t) + m\{2m\ - s - u) + m~u - tu))))) 

/((i-m|o)(i-m|o)) 

|r|2 = ^TJxTJ (B20) 

Channel II (the i-channel I exchange) is now summed over Ii and I2, with appropriate 
propagators. Also added is the neutralino w-channel exchange. 

= -g2/V2{cos9e{-Ni2-tein9yyNii)+sm9eme/{myyCosP)Ni3) 

fAi2,i) = -g2/V2{- sm9i{-Ni2 - tanOwNii) +cos9eme/{mw cos P)Nis) 

fsC^^i) = -g2/V2{cos9im(,/{mw cos P)Ni3 + sm9i{2tan9^vNii)) 

fB{2,i) = - g2/V2{- sin 9ime/{mw cos P)Ni3 + cos 9i{2 tan 9w Nil)) 

= -Wcose^((-l/2)cos2e,-(-l)sin2M 

/c(2) = (72/ cos ^^((-1/2) cos sin ^^,) 
h = -Wcos^^,,.((-l/2)-(-l)sin2^^) 
fn = ^2/cos6'vv((-l)sin2^,^) 
/olW = g2/i2cos9^)iNi^Nu- NMs) 

foR{i) = -g2/{2cos9^^){N,^Nii-Ni3Ni3) 

TixTi = {l/2){2fA{l,l)fB{l,l)m^me{2mtfl + QmjmlfLfR + 2mtfl 

+ml{fl + fn) {m\ - .) + m\ {fl + fl) {m\ - s) - 2ml f Is + Qmlhhs 
-2ml fls - -mjflt - myit + fist + f^st - mjflu - m]fl.u + flsu + f^su) 
+/a(1, lf{-m]m\fl + m\jl{m] - s) - m\fls - QmjmlfLfRS 
+m^m|/^s + mlfls'^ + mpK-m'j + s) - m'jmlf^t + m'jflst 
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I 2 t2 A. £2 24. , A £2 I 2 2 ^^2 2 ^^2 2 ^^2 

+m|/i(2mt/i + mliGmjfR - 2/^^) - /^K' -s)is + t + u)) 
+m\{~2m%fl + 2mlfL{-^m]fR + f^s) - {m\ - s){m)Jl - fl{t + u)))) 

l)^(mj/^(m^ - s) - mjflu + m]{mifl - m\jl. - f^st + flsu 
+'^\Jr{s + t + u)- mlfL{^m]jR - 6fRS + /^(s - t + u))) 
+m\{m\fl + m]{Qmlf,fR - fh - flit + u)) + fl{2ml 
—2rnls + s{t + u))) — /^(— m|^s — s(m^ + st — m\{s -\-t — u)) 
+m|(2m^ - 2mls + s(s + i + u)))))/{ml{s - mff) 
TuyTu = (l/2)((/^(?,l)(-2/B(j,l)m^m, + /A(j,l)(m| + m2-t)) 
l)(-2/A(j, l)m^me 
+fB{j, l)(m| + mj - t))){ml + - 2mlt + - 2m\{ml + t))) 
/(ml(t-ml)(t-ml)) 
^iixT^n = l/2((l/m^)(/B(l,i)(/^(l, j)m£(m|(m^o/oL(«)/oL(j) + m^ofoR{i)foR{j)) 
+m\{m^ofoL{i)foL{j) + m^ofoR{i)foR{j)){ml - 2u) 
-6rn^ml{m^om^ofoR{j)foR{i) + foL{i)foR{j)u) - {m^ofoL{i)foL{j) 
+m^ofoR{i)foR{j)){'2rnl - '^l'^ - u^)) + fB0-,3){m\foL{i)foL{3) 
-4:m]mlfoL{.'i)foL{.3) - '2mlfoLii}foLij) - rn^om^omlfoR{i)foR{.3) 
-3mlml{m^ofoL{j)foR{i) + m^ofoL{i)foR{j)) + '2mlfoL{i)foL{j)s 
+m^om^omlfoR{i)foR{3)s - ^rn^rnl{m^ofoL{j)foR{i) 
+m^ofoLii)foRi3))i'^m'^e + ml - s - t) + 2mlfoL{i)foL{,3)t 
-m^om^omlfoRii)foRi3)t + "^J/oL(0/oL(j)(2"^| + 2m| - s-t - 2u) 
+2rr4mlfoL{i)foL{j)u + 2mlfoL{i)foL{j)u + rr4m^om^qfoR{i)foR{j)u 
+m^om^omlfoR{i)foR{3)u - '^rnlfoL{i)foL{j)su - m^om^ofoR{i)foR{j)su 
+myoL{i)foL{3)u^ - foLii)foL{j)tu'^ + rn\{-mlfoL{i)foL{,3) 
+m^qm^ofoR{i)foR{j)s + foL{'i)foL{j)su + 2foL{i)foLi3)tu 

1 3 

+foL{i)foL{jW - mlfoL{i)foL{j){s + t + u)+ rr4{2mlfoL{i)foL{j) 
-m^om^ofoR{i)foR{j) - ^foL{i)foL{3)u)))) 

* 3 



59 



i){fB{'i-,j)mi{m^{m^ofoL{i)foL{j) + m^ofoR{i)foR{j)) 
+ml{m^ofoLii)foL{j) + m^o/oR(i)/oK(j))(m| - 2u) 
-Qm^ml{m^om^ofoL{i)foR{j) + foL{j)foR{i)u) - {m^ofoL{i)foL{j) 
+m^ofoR{i)foR{j)){'2rnl - '^l'^ - w^)) + fAO-,j){m\foR{i)foR{j) 
-Smlml{m^ofoL{j)foR{i)+rn^ofoL{i)foR{j)) - ^rn^rnl{m^ofoL{j)foR{i) 
+m^ofoLii)foRij))i'^m'^e + ml - s - t) 
+myoR{i)foR{j){2mj + 2ml-s-t- 2u) 

+m^om^ofoLii)foLij){-ml + ("^1 " ^)'^ + "^li^ -t^u)) 
+foR{^foRij)i-2ml - 2mlsu - tu" + 2m\{s + t + m) 
+m'j{-Aml + 2mlu + i?)) - m'^~{m,lfoR{i)foR{.3) 
-m^om^qfoL{i)foL{j)s - foR{i)foR{j)su - 2foR{i)foR{j)tu 
-foR{i)foR{jW + 'mlfoR{i)foRU){s + t + u)+ mj{m^om^ofoL{i)foL{j) 
+foR{i)foRm-2ml + 3u)))m/{{u - m|o)(« - m|o)) 
(l/2)(-/^(l, 1)(2/b(j, l)m^me{mlh + 2mlh - SmjmlfR - m\fR 
-2ml fLS + fnlfRS + rn^fLimj^ + 3ml - t) - Sm^/it + mjfRt + mlfRt 
-fRst + ht'^ - mlfiu + 2ml /ru + htu + m|(m|(/i, - 3/r) - m|/R + /rs 
-'2 fit - !lu)) + /^(j, l){-2m\m\fR + mfmlfL-s + mlfis + AmjmlfRS 
-mlfLS^ + mjmlfLt + m\fLt + 2m]mlfRt - m^fLst - mlfLst + JlsH 
-mjht'' - 2ml fLt" + fLt^ + mlfL{-2mj + s + t-u)- mjml^u 
-m^fLU + 2m'jmlfRU + m'jfitu + mlfitu - 2m'jfRtu + mlfiu^ - fitu^ 
+m|(m|(2m|(/i - 3/,^) + his + 3t - u) + 2fRu) + h{-s'' - st - 2^ 
+2ml{s -t-u) +tu + u^)) + m\{2ml{fL - /ii)(3m| - t) 
+m|(2m|(/L - Jr) + /i(8m| -s + t-u))- hitis + t-u) 
+ml{s - t + 5m))))) + l)(-2/A(j, l)m^me{-mlfL + mjjR + 3m\mlfR 
+2mlfR + mlfis - 2ml fns + ^IfLt - mlfRt - 3ml JrI - fhst + /^t^ 
+m'jfL{-3ml + t) + 2ml Jlu - '^Ifn'^ + fnt'^ + {-"^Ih + m\fR 
+mli-3fL + /r) + hs - 2fRt - fRu)) + /bU, l){2mWzfL - ^m^mlhs 
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-u^mlfRS - mlfRS + mlfRS^ - 2rrvlm\fLt - m^mlfRt - mlfRt 
+m'jfRst + mlfRst - JrsH + m]fRt^ + 2m\fRt^ - fRt^ - 2m'jmlfLU 
+m]mlfRU + m\fRU + 2m]fLtu - mjfRtu - mlfRtu - mlfRU^ + fRtu^ 
+mjjR(2mj -s-t + u) + m\{-fR{-s^ - st - 2t^ + 2m|(s -t-u)+tu 
W) + mj{ml{6fL - 2fR) - 2hu + - 3i + u))) + m|(2m|(/i 

-/K)(3m| - t) + m\{2m]{fL - Jr) + /i?(-8m| + s-t + u)) + fR{t{-s 
+i - m) + ml{s - t + 5u))))))/{2ml{s - ml){t - m|)) 
= l/2((l/(2m2j)(/A(l, 0(2/^(1, l)m,(m|/R/oR(i)+m|m^o/oi(i)(-m|/L 
-^mlfR + /^s) + m|/oR(2)(m2 (-3/l + Jr) + /,j(m| - t - 2u)) 
+rn^foR{i){-Q'rnlmlfL - ml{3fL + ^Jr) + ml(-3mljR + 3/l(s + t) 

+ 2t)) + fRu{-m^^^ + t + u)) + m^ofoL{i){-mlfL + - sjii 

+m^(3m|/R - 3/rs + /^(s - t + li)))) + /^(l, l)(2m|m^o/i/oi(i)(-m| + s) 
+m\foR{i){2mlfR - fii^s + t-u))- 2m^m^ofoLii)i2mlfL + m}jL{rnl - s) 
+ml{6mjfR - 2hs) - himj -s){t + u)) + m|/oR(z)(m|/L(4m,2 + 8ml 
-3s-t + u)+ ml{2mlfR - AfRU + - 3t + u)) + /^(s^ + f + tu 

-2u^ + s(2t + u)- m|(5s + t + 3u))) + foR{i}{mjjL{4^ml - ^rn\{s + t) 
^m\{8m\ - 2m) + (s + t - u)u) + m\{-Am\fR + ml{-8fLS + 2/rm) 
+u{fL{s + t-u) + 2Jru)) + /L(-4mts - uis" + e - u^) 
+ml{As'' + {t- u)u + s{At + u)))))) + i)(2/A(l, l)m,{m\fLfoL{i) 
+m\m^ofoR{i){-3mlfL + fRi-rrij + s)) + m^foLii){.-m%{AfL + 3//?) 
+m|(-6m^/H + 3/ij(s + t) + /l(s + 2t)) - m|/L(3m^ + + /^^(t + u)) 
+rn\foL{i){m]jL + - 3/i?) - fiit + 2u)) + rn^ofoR{^){^rnj^rnlfL 

-mlfR + fR{ml - s)u + m|(-3/LS + - t + i*)))) 

l)(2m|m^o/jj/ofi(i)(-m^' + s) + m|/oL(i)(2m^'/L - + t - «)) 
-2m^m^ofoR{i){mj{6mlfL + - t - ti)) + fR{2ml - 2mls + s(-m| 

+i + u))) + mlfoL{i){mj{2mlfL + Am]jR - /rs - 3fRt - Af^u + /ru) 

+ 2si + + su + tu- 2u^ + m| (8m^ - 3s - i + li) - m|(5s + 1 
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+^u))) + foL{i){m1{-AmifL + 2m|(4m|/R - Afns + hu) + u{-2mlfR 
+fB.{s + t-u) + 2fLu)) + /R(-4mts + m| (4mt - 4m|(s + t) + (s + i 

- ^(s^ + - u^) + m|(4s^ + - u)u + s(4t + «)))))))) 
/ {{u - m\o){s - m])) 
TJixTJii = l/2((l/m|)(/B(j, 1)(/a(1, i)mi{m^ofoL{i) + m^foR{i)){m\ - 2m\s 
+m'^{mj^ + 3m| — t) — + (3m| — -u) — m|M + tu) 

i)(mJ/oL(i)(m| + 3m^ - ^) + m^m^ofoR{i){-ml + mis 
+m|^(— — 3m| + s) + m^t — st + 'm'j{—3m1 + t) + 2'mlu) 
+rn\foL{i){'iml - Amis + m^(6m| - 2t) - 5m^t + si + + m|(2m| 
+5m^ — s — t — m) — 2m^u + tw) + foiivi—'mz-^ + ^z-^ ^ ""^z^ + 2m2st 
+mlt^ + m\u + m|iw — t^u — mlu^ + m^(— m^(3s + + iu) 
+m|^(m|(4m| — li) + + ml{—s — 3t + u))))) 

l)(/ij(l,i)m£(m^/oL(«) + m^^^foRii))iml - 2m|s + m|(m| 
+3ml — t) — mlt + mj^{3ml — u) — mlu + + /a(1, i)ifTi^foRii)i^i^ 
+3ml — t)+ m-^m^o foL{i){—fn\ + "^^s + ~ 3m| + s) + 

—si + m'^{—3ml + t)+ 2mlu) + m|/oR(i)(3m^ — 4m|s + mf{6ml — 2t) 
—bmlt + st + + (2?Ti£ + 5m| — s — t — m) — 2mlu + tw) 
+/oj?(i)(— m^s + m^s^ — m^t + 2mlst + mlt^ + m^u + m|iM — t^u — mlu^ 
+m'j{-ml{3s + t)+ tu) + m|(m|(4m^ -u)+tu + m^(-s - 3t + u))))))) 
/((tz-m|o)(t-m|,)) 
|T|2 = %xT, + J2fc{^fc{j)TuxTu + J2^uixTui + 2j2fc{mxTn 

+2 ^ TJxTin + 2 ^ fc{j)TuxTjn (B21) 

^iX — ^ 7^ 

The couplings are modified. 

/a = -g2/V2{cos9i{-Ni2-tein9wNu)+sm9ime/{mwCosP)Ni3){e) 
fs = -5f2/>/2(cos6l£m^/(mtvCOS/3)A^i3 + sin 61^(2 tan 6'w^7Vii))(e) 

62 



TiixTu = {l/2){-2{-4:fAfBm^mi + fl{ml + mj-t)+fUrnl + mj-t)){ml+t)) 
I it -ml? 

TixTii = (l/2)(l/2(-4/A/ijm^m^(3m|-3m| -3m2-s-t + M) + /2(-5m2s + s2 
—rnjt + t'^ — mju — + {2ml — s + 3t + u) + m|(— 8m|^ + s — i + 5m)) 
+/!(— 5m^s + — mft + t'^ — rnju — + (2m| — s + 3i + w) 
+m|(-8m| +s-t + 5u))))/{{s - n4){t - m|)) 
|T|2 = 71x71 + TIixTIi + 27Ix7Ii (B22) 



The s-channel r annihilation was neglected in [17] due to the small Yukawa coupling. 
However, at large tan /3, these couphngs are enhanced particularly for the h final state. The 
t-channel ^exchange is now summed over Ii and £2- Also added is the neutralino li-channel 
exchange. 

I. s-channel i annihilation 

II. t-channcl £1^2 exchange 

III. ti-channel Xi,2,3,4 exchange 

/i = —g2iTit cos a I {2mw cos jS) 



—milt — su — m^{mi +t + u)+ m|^(2m£ + s + t + u)) + Jev^x ~ '^h 

— 3m^s — m^t — su — m^{rn1 + t + u) + mj^{2m'j + s + t + u))))/{s — m^Y 



/A(l,i) 



/2(1) 



/2(2) 



g2mz{{—l/2 + sin^ Ow) sin[— cos](/9 + a) cos^ Og, 

— sin^ ^vv' sin[— cos] (/3 + a) sin^ 9() / cos ^v^. 
-\-g2ml sin[— cos]a/(mvi. cos /5) 

— fi'2^£(^£ sin[— cos]a — /x cosfsinja) sin 6'^ cos 6'£/ (mi^ cos /3) 
92mz{{l/2 — sin^ ^^y) sin[— cos](/9 + a) 

— sin^ 9yy sin[— cos](/3 + a)) cos ^£ sin 9e/ cos 

— 5'2'^£(^^ sin[— cosjo; — /i cosfsinja) cos(2^^£) / (2m;v cos /9) 
-g2/\/2{cos6e{-Ni2 - tanOwNn) + sm6ime/{mw cos/3)7Vj3) 
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/A(2,i) = -g2/V2{-sm9e{-Ni2-tein9wNii)+cos9emi/{mw COS P)Nis) 
fB{l,i) = -g2/V2{cos 9ime/{m^^. cos (3)Ni3 + sm9e(2taii9^Nii)) 
fB{'2,i) = -g2/ V2{- sm9eme/{mw cos P)Nis + cos9i{2tan9wNa)) 
fcL{i) = y2/2((A^i3(A^i2 - A^iitan^;^) + Ni^{N^2 - A^iitan^^^)) sin a 

+(iVi4(A^i2 - A'"*! tan 9^) + Nii{Ni2 - Nu tan 9^)) cos a) 
fcRii) = ^72/2((iV,3(iVi2 - iVntan^^^) + iVi3(iV,2 - iVatan^,^,)) sin« 

+{Ni4{Ni2 - Nu tan 9^^) + Nu{Ni2 - Na tan 9^)) cos a) 
%x% = l/2(/^(l,l)(2/B(l,l)m^m^(m| + m| + 3m^' + s-i-ii) 

— + + 3m^s — m^t + m'^{5m'^ — t — u) — su 
+ml{-Amj + s + t + u))) + 1)(2/a(1, l)m^m^(m| + m| + Sm^ + s 
— i — + — + + 3m^s — m^t + m|(5m^ — t — u) — su 

+m\ (-4m^ + s + i + - ml)^ 

TiixTii = l/2(/^(^,l)(2/B(jU)m;^m, + /^(j,l)(m| + m,^-t)) 

1)(2/^(j, l)m^m, + l)(m| + - - m|)(t - m|,)) 

^iixTJii = 1/2(/a(1, j)rne{fcR{i)rn^o{2fcL{j)rn^m^o + fcR{j){m\ - ml^jj^ + u)) 

+fcL{i){'2fcR{j)m^u + fcL{j)m^o{ml - m^H] + 

j)(/ci?(^)"^x?(/ci?(j)"^x;'("^l + m]-t) + fcL{j)m^{m\ + m^[^] 
+2m^ -s-t)) + fcL{i){fcR{3)rn^rn^o{m\ + m^[^] + 2m^ - s - 

- ^Mi/] + ^|(2m| - s - + (m| - + m^[^](-2m| + s + i 
+«))))) + fB{l,i){fA{l,j)me{fcL{i)m^o{'2fcR{j)m^m^o + fcL{j){m\ - m^H] 
+u)) + fcR{i){2fcL{j)m^u + fcB.{j)mxo{ml - m^jj^ + 

+fB{l,j){fcL{i)m^o{fcL{j)m^o{ml + mj - t) + fcR{3)m^{m\ + m^[^] + 2m| 
-s - t)) + fcR(i)(fcL(j)m^m^o(ml + m^H] + - s - + fcR{j){m\ 
-iTT'hiH] + m\{2m\ - s-t) + {m\- s)u + m\ji^{-2m\ + s + i + «)))))) 
/((-u - ?n|o)(M - ?«|o)) 
T^xTJi = 1/2(/b(1, l)mt{2m\ - m\ ^m\^s-t)^ fA(i, l)m^(m| + m| 

+Amj -t-u)) + l)m^(2m| - m| + m| + s - 
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l)m^{ml + m\ + Am] - t - u)))/{{s - rr4){t - m|)) 
TixTni = 1/2((1/2)(/a(1,2)(2/a(1, l)m,(/cR«m^o(2m| - m| + m] + s - t) 

+fcL{i)m^{ml - m\H] - 2?n| + 2m^ + u)) + fsil, l)(2/ci?(i)m^m^o(m| 
+m|^ + Am] — t — u) -\- fcL{i){'2m'^ + mj^s + m]s — — rnj^t — m]t + 
+m|(2m^[^j + 2m|^ + Gm] — s — 3t — u) + rrij^u + Smju — + m^[^](— 2m|^ 
-2m] + s-t + u)))) + 0(2/^(1, l)me{fcL{^m^o{2ml - m| + m] + s 
-t) + fcR{^)m^{ml - m^H] - 2m| + 2m] + u)) 

l){2fcL{i)rn^rn^o{rn\ + m| + Am] - t - u) + fcR{'i){'^rn\ + m|s 
+m]s — — m|^i — m^i + t'^ + m^^{2m\^-^^ + 2m|^ + Qm] — s — 3t — u) 
+m'j^u + 3m^M — + m1^ff-^{—2mj^ — 2m] + s — t + «)))))) 
/((s-m^)(K-m|o)) 
TJixTJii = 1/2(/b(j, i)(/cL(0^x°("^l + ^1 - ^) + fcR{i)m^{rn\ + m^^^j 

+2m^ - s - t)) + fA{l,i)mi{2fcR{i)m^m^o + fcL{'i){m\ - mly^^ + m))) 

i){fcR{i)m^o{ml + m] - t) + fcL{i)mx{m\ + m^[^] + 2m] 
-s - t)) + i)me,{2fcL{i)m^m^o + fcR{i){m\ - m\jj^ + i*)))) 
/((t-m|,)(^-m|o)) 
|r|2 = /27Ix7i + X:/2(^)/2(j)^ix7Ii + 5]7Inx7In + 2^/i/2(0Mi 

+2 ^ /iTJxTin + 2^/2 (jO^iX'T^ii (B23) 

I. s-channel £ annihilation 

II. i-channel I2 exchange 

III. w-channel Xi,2,3,4 exchange 

= -g2/V2{cos9e{-Ni2 - t&ney^^Nn) +smeeme/ {m^ cos P)Ni3) 

/^(2, z) = -5(2/ V2(- sin6'£(-iVj2 - tan^^vv-iVa) + cos6emi/{mw cos P)Ni3) 

= -g2/V2{cos6emi/{mw cos P)Ni3 + sm6i{2tan6wNii)) 

fB{'2,i) = -g'2/v^(-sin6'^m^/(mvvCOS/3)A^i3 + cos 61^(2 tan 6';vA^ji)) 
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fi = -{g2'meteiiiP)/{2mw) 

fciii) = g2/2{{Ni:i{Nu-Nnta.ne^) + Nis{Ni2-Natane^))sm(3 

-{Ni4{Ni2 - Nn tan0^) + Nu{Ni2 - Na tan^^^)) cos/3) 
fcnii) = -92/'2{{Nis{Nr2 - Nn tan 6^) + Nis{Ni2 - tan 9^)) sin /3 
-{NiiiNi2 - Nn tane,,,-) + Nu{N2 - Ni tan 9^)) cos /J) 
TixTi = l/2{fA{l,){2fB{lA)rn^m{-rnl-ml-mi + s + t + u) 

1)(— m~ + — mj + m^s + mjt + su + nt^{—mj + t + u) 
-m|(s + t + + 1)(2/a(1, l)m^m^(-m| - m| - + s + t + li) 

1)(— m| + — + m^s + m^t + su + m'^{—ml + t + u) 
-ml{s + t + u))))/{s-mj)^ 
TiixTu = l/2(/A(2,l)(2/B(2,l)m^m^ + /A(2,l)(m| + m|-i)) 

+/b(2, 1)(2/^(2, l)m^m, + /b(2, l)(m| + m,^ - t)))/{{t - m\f) 
= 1/2(/a(1, i){fB{lJ)me{fcR{i)m^o{2fcLij)m^m^o + fcR{j){-ml + m| + «)) 
+fcL{i){'2fcR{j)m^u + fcL{j)m^o{-ml + m| + w))) 

+fA{l,j){fcR(i)m^o{fcR{j)m^o{ml + mj-t) + fcL{j)m^{m\ + m| + 2m| 
--s - t)) + fcL{i)ifcRij)m<^m^«iml + m| + 2m^ - s - t) + fcLij){-m\ 
+mj + m\{2ml - s - t) + {mj - s)u + m\{-2ml + s + i + «))))) 

i){fA{iJ)rne{fcL{i)rn^o{2fcR{j)rn^m^o + fcL{j){-ml + m| + w)) 

+fcRii)i'2fcLij)mxU + fcRij)m^o{-ml + m| + m))) 

+fB{^J){fcL{i)m^o{fcL{j)m^o{ml + mj - t) + fcR{j)m^{m\ + ^ + 2m^ 
-s - i)) + fcR{i){fcL{j)m^m^o{m\ + m| + 2m| - s - + fcR{j){-m\ 
+m| + m|(2m| - s - + (m^ - + m^(-2m| + s + i + w)))))) 
/((« - m|o)(M - m|o)) 
TJxTJi = 1/2(/a(1, 1)(/a(2, l)m^(-m| - m^' + s + t) + /b(2, l)m^(m| + m| - t - «)) 
1)(/b(2, l)m^(m| + mj - s - t) + /a(2, l)m;f(-m| - m| + t + u))) 
/{{s-mf){t-ml)) 
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TixTJii = l/2((l/2)(/B(l,i)(2/B(l,l)m,(/cL(iKo(m|+m|-s-i) 
+fcR{i)m^{3ml + m| + 2m| + 2m^ - 2s - 2t - «)) 

l)(2/cL(«)m^m^()(m| + m| - t - m) + fcR{i){'2rn'^ + m\s + m^s 
— — — m^i -\-t^ -\- m^{2rnq^ + 2m£ — s — 3i — li) + rin^i^u — ntfu — 
+ml{2ml - 2m| + 2m| + s - i + -u)))) 

+/a(1,0(-2/a(1, l)"^£(/ci?(«)m^o(m| -s-t) + fcLii)m^i3ml 

+m| + 2m| + 2m^ - 2s - 2i - «)) + /^(l, 1){2 fcR{i)rn^m^o{-m\ - m| 
+i + + /cL(0(~2m| — m|^s — m^s + s^ + m|^i + — — rnj^u + mju 

- m\{2m\ - 2m| + 2m] + s-t + u)+ m|(-2m| - 2m] + s + 3i 
+u))m/{{s-m]){u-m^.)) 
TuxTui = 1/2(/b(2, i)(/cL(i)m^o(m| + - + /cii(i)m;^(mi + m| + 2m] 

-s - t)) + i)me{2fcR{i)m^m^o + fcL{i){-m\ + m| + u))) 

+/a(2, l)(/A(l,i)(/ci?,(^)"^x°("^l + - t) + fcL{i)my.{m\ + m\ + 2m] 
-s - t)) + fB{y,i)me,{2fcL{i)my^m^o + fcR{i){-m\ + m| + «)))) 
/((i-m|)(ii-m|o)) 
|T|2 = f^TjX% + f^TuxTn + J2TjnxTni + 2f,f2TixTn + 2j2fi'rixTni 

hj i 

+2^/2^1X^11 (B24) 

Appendix C: Chargino-Slepton Coannihilation 

Below is the hst of the amphtudes squared for chargino-slepton coannihilation. Note that, 
for identical-particle final states, one needs to divide them by two when performing the 
momentum integrations. 

hxi — ' ^cZ 

I. i-channel xji 2) exchange 

II. u-channel ^(l^2) exchange 
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flLR{i) 


/^L /^R 


fiRhii) 


/^R /^L 

- ^~i^-vt-x- ^xt-xt-z 


flRR{i) 


r<R nR 


f2L{i) 


— _ 


f2R(i) 






= {'m-+{3fiLL{j){flLR{ 



+m^+{fiRL{i)fiRL{j){ml + s{-ml+ +t)- ml{s + t-u)) 
+fiLRii)fiLRij)iml + s{-m'i+ +t)- ml{s + t - u)))) 



+/iLL(0(3(/iLii(j) + fiRL{j))m-+m-+ml{mj^ - t) 

n / \ / A A O 00 O 

— 2/ii,l(j j \—m^t — m^t — m^st + m^t + m^tu — t u 

+m| (2m^ - 2mlt + st + t'^ + tu) + m|+ (m| + - m| (s + i + h))))) 

/(m|(mL-i)(-mL+i)) 

-2mf {ml + u)))/{ml{-ml + u){-ml + u)) 



^1 

TJxT^i = i-2{fiLRij)f2Lii) + fiRL{j)f2R{i})m~+m~+{ml + m|(3m| - s) + s-u 
-m^(s + 2i + ^)) + fiLL{j){f2L{i) + /2rW)K\(s - t - li) 



+ (m| - u){-s^ + + ml{s -t-u)+u^) 
+m|(-s^ + - su + tu + 2u^ + 2771^(5 -t + u)) 
+m?+(2m|^ - (m| - u){-s + i + - m|(2m^ - s + i + 3ii)))) 
/(2m2(m2+-t)(m|-tz)) 
|r|2 = ^(7ix7i + 7iix7ii + 27ix7ii) (CI) 



^1X1 — ^ i^^7 

I. t-channel Xi exchange 

II. ti-channel £1 exchange 



•^1^ - ^ii-..-xr ^xt-xt-i 

fiR = C~+_j+_^ 
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f = c~ ~ ~ 

TlXTi = (2(/i\ + (mj+ - m| - st - m|+ (6m| - 5t + + m| (s + t + «))) 



/(^+ - tf 



TnxTu = -2(/2\ + /|K)(m|+-«)(m|+^/)/(m|-^x)' 
^xTii = {{f,Lf^^ + f,^f^j,){-s^ + t^ + u^-ml+{6ml-3s-3t + u) 
+m| (s - i + 3M)))/(2(m|+ - i)(m| - u)) 
|r|' = TixTi + TuxTu + 2TixTu (C2) 

I. t-channel Xi,2 exchange 

II. M-channel ti,2 exchange 

hLR{i) = cl_^^_~.c^+_~+_^ [c|_^^_~_ c|+_j+_^] 



TixTi = {m-+{fiLL{i}ifiLR{j)rn-+{-ml +t) + /iLL(i)m-+(-m?+ + u)) 
+fiRR{i){fiRL{j)m^+{-m\ +t) + fiRR{j)m^+{-ml+ + u))) 
+fiLR{i){fiLL{j)m^+m-+{-mj^ +t) + fiLR{j){mj^ + m|+(2m| - t) 
+st -mj^{s + t + u))) + fiRL{i){fiRR{j)m^+m^+{-mj^ + t) 
+fiRLij)iml + m\+{2m\ -t) + st- m|(s + t + u)))) 

^ixTii = (/2^(z)/2^(j) + /2«(^)/2K(j))(m|+-«)/(-m|+«)(-m|+«) 
= {fiLR{i)f2L{j)m^+{mj^ - t) + fiRL{i)f2R{j)m^+{mj^ - t) 

+(fiLL(i)f2L(j) + fiRR{i)f2R(j))m-^{ml^ - u))/((ml^ - i)(m| - u)) 
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|r|2 = ^(7Jx7i + 7Jix7Ji + 27Jx7ii) 



I. i-channel Xi,2 exchange 

II. li-channel £2 exchange 

fiLLii) = Ci_,^_~_Ci+_~+_^ 
flRRii) = Cf+_~+_^ 

/2L = C'^,.^.?'- 

f2R = Cj^.A-l-C~--.,-J, 

TixTi = {m^+{fiLL{i){fiLR{j)m^+{-mj^ +t) + fiLL{j)m^+{-ml+ + m)) 
+fiRR(i){fiRL{j)m^+{-mj^ +t) + fiRR(j)m-+{-ml+ + u))) 
+fiLR(i)(fiLL{j)m^+m^+{-mj^ +t) + fiLR(j)(mj^ + m?+(2m| 
+st - m|(s + t + u))) + fiRLii){fiRRU)m^+m^+{-mj^ + t) 
+fiRL{j){ml + ml+{2m\ - t) + st - m\{s + t + u)))) 
/{{ml.-t){-ml,+t)) 
T,,xT,, = {fL + f2R)Ht-n)/H,-^r 

TlxTii = {fiLRii)f2Lm^+imj^ - t) + /iKL(i)/2/?m-+(m| - t) + (/ill(0/2- 
+fiRR{i)f2R)m-+{ml+ - u))/{{ml+ - t){m^^ - u)) 



I. s-channel annihilation 

II. i-channel 2,3,4) exchange 

Note the L-R switch for the neutrahno couphngs below. 
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flR = ^i-^.,-^- Cue-e-w+ 

f2LL{i) = Cf^_i_^o C^_^-_^o 
f2LR{i) = C'|_^_~o C'^_J-_J0 

f2RL{i) = Cf^_e^~oC^_~-_~o 

f2RR{i) = C^^_^_~,C^ 

TjxTi = 1/ (m^s^) (/i^,(m|_^ {mj - s) + mt+ (-m^ + s) + s(-mt + (m| - s)u 
+m^(s - i + m)) + m^+(-2m^ + 2m^s + (m^ - + u)) 
+m|(2m^ - 2m^s - (m| - + i + w)))) 
^iXT^i = irn^o{f2RLii)i3f2RRij)rn^+ml,{mj^-mj-t) 

+m^o{6f2LRij)m^+miml. + /2RL(j)("^t + - s) + + s)t 

-m^(s + i - u)))) + f2LL{j){^f2LR{i)m-+ml,{mj^ - - t) 
+f2RL{i)mi{-2rn^ + 2m^i + i(m| + - s - + rni+{-m'~_^ - ml 
+s + m))) + f2LRii)im^oiQf2RLij)m^+meml, + /2Li?(j)("^t + ^l+i^'e - s) 
+ {-mj + s)t - m^(s + t - m))) + f2RR{j)mi{-2m'^ + 2m^t + t(m| + m| 
-s - li) + ml+{-m'~_^ - + s + i*)))) + f2LL{i){^m^+mlr{f2LRij)rn^o{mq^ 
-m] -t)+ 2f2RR{j)met) + f2RL{j)m'm^o{-'2'rn% + 2m^t + t(m| + m| 
-s - m) + ml+{-mj^ -ml + s + u)) + /2LL(j)(2m^m^ + - m'^t 
—2rrilm\^t + m^i + mist + m^^st — m^t^ + m^tu — m^tu + t^u 



+to|+(— + — tu — m|(s + + + i + li)) — mj^{2m:^ 



4 

2 



-2m^t + t{s + t + u)))) + f2RR{i){3rn^+m^.{f2RL{j)m^o{m^^ -m^^-t) 
+'^f2LLij)mt) + f2LRij)mi'm^o{-'2ml, + 2m^t + t(m| + - s - m) 
+m^+(-m| - + s + + /2iiii(j)(2m|m^ + m|^i - m^i - 2m|m^t 
+m^f + m|si + m'^st — m^t^ + mjtu — m'^tu + t^u 
+m|+(— + — tu — m'j{s + u) + Tnf^{s + t + u)) 
-ml {2m% - 2m\t + t(s + t + u))))) 
/{ml{mlo-t){-mlo+t)) 
TixTii = ifiRi-2f2LL{j)m~+mi{mj^ + 4m^ - m^s - - m^i + si 
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+m|^ {rrij + 4m^ — s — t — 2u) + m|+ (m|^ + mj — s — u) — mju 
-Sm'i.u + su + tu + u^) + 2mj^o(3/2LR(j)m£m^(m|+ - m| + s) 
+f2RL{j)m-+ (2m^ + m|+ (m| - s) + m| (m^ - s) - 2m\,s - mjt + st 
-mju + su)) + f2RR{j){-5m'jml,s - m^s + m^s^ - m^ml^t - m^t 

n 2 4. , 2 4.2 4/ij. \| 22 , 4 , 2 , 2 

—Irriy^st + m^i — ''77.^^(5 + i — it) + m^rriy^u + m^it + m^srt + rriy^su 
—s^u + m^t-u + rn^tu — t^u — m^-u^ — 2m\Ai^ + + m?+ ((— 2m£ + s 
+i — u)u + m^(— s — 5i + m) + (4m£ + 8m^ — 3s — i + «)) 
+m| (4m^ + + 2si + + + - 2u^ + m]{d>ml^ - s - 3i + 
-2m^(s + i + «)))))/(2m^s(m|o - t)) 
|r|2 = 7Ix7i + ^7Iix7Ii + 2^7Ix7Ii 



I. s-channel annihilation 

II. t-channel 2,3,4) exchange 

Note the L-R switch for the neutralino couplings below. 



/2Li?(i) = C'|_^_~,, C^+_~-_~o 



/2i?L(i) - C~_^_~o C^+_~-_J0 

TixTJ = (/2^(^^_^|+s(m,'-t)+m^+(2m|-i-ii) 
+m|(-2m^ + s + t + m)))/s2 
%ix%i = (m^o {f2RR{i) {f2RL{j)m^+{-mj^ +m'j + t)- m^o {2f2LLU)m^+me 

+f2RR{j)iml+ +mj- u)) + f2LR{j)me{2m'i+ +mj^+mj-s- u)) 



+/2LL(0(/2Li?(j)"^^+(-"^| +m'j + t) - m^o{2f2RR{j)m^+me 
+/2LL(i)(m^+ +mj- u)) + f2RL{j)me{2m'i+ + m| + - s - u))) 
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+f2LR(i)(f2LL(j)m-+m^Q{-m^^ + + t) + me(-2f2RL(j)m~+t 
+f2RRij)m^o{2ml+ + mj^ + nij - s - u)) + f2LRij)imj^ - mj + mjs 
+m|+ (2m|^ — 2rnj — t) + st + mju — (s + t + u))) 
+f2RL{i){f2RR{j)m^+m^o{-mj^ + + t) + me{-2f2LR{j)m^+t 
+f2LLij)m^o{2m'i+ + m| + - s - u)) + f2RLij)imj^ - + mjs 
+m|+ (2m| - 2mj - t) + st + mju - m| (s + t + «)))) 
/((m|o-i)(-m|o + i)) 
TixTji = -{fiR{m^o{f2RR{j)mi{m'i+ - m| + s) + f2LL{j)m~+{rni+ + m| 
+2m| - i - - f2LR{j){mj^ + m^^+{mj^ + mj) + st 
+mj^{mj - s - t-u)) + f2RL{j)m-+mi{-2ml+ - 2mj + s + t + u))) 
/{s{mlo-t)) 

|r|2 = TixTj + ^rnxTn + 2^TixTn (C6) 



Iixt w-e+ 

I. t-channel z?^ exchange 

fiR = C'i!.j;^„^- C~ 
TixTi = ii-AfiLfmm^+me + fiL{'m\+ + - t) + flR{m\+ + m] - t)){m\ 
+ {ml - tf - 2ml{ml + t))) / {ml{m% - tf) 
|r|' = TixTi (C7) 

I. t-channel exchange 

fiL = C^_~^_~- C~_7^_^+ 

%xTi = (-4/ii/,^m-+m, + A\(m^++m,2-t) + A2«(m^++m,2-i))/(i-m?j2 
|r|2 = TJxTi (C8) 



73 



4xr w- 

I. t- channel 234) exchange 

II. K-channel exchange 

Note the L-R switch for the couphngs below. 

fiLR{i) = C~_e_^oC^-^-_~o 

flRL{i) = Cj ~o 



flRR{i) — Cj ~o 



f /^R ^ ^ „ 

f2R = ^i-Vt-xT ^vt-h-W 

TixTi = (m^o(/iRL(i)(3/iRi?(j)"^-+m^(m| - - t) 

+m^o(6/iLR(j)m-+m^m^ + /iRL(j)("^t + ^l+i^I - s) 
+{-ml + s)t - m^(s + t - u)))) + fiLL{j){^fiLR{i)m^+ml,{mj^ - 
+fmL{'i)me,{-2rn^^ + 2m^t + i(m| + - s - li) 
+m^+(-m| - + s + m))) + fiLR{i)im^o{6fiRLU)m^+meml, 
+fiLR{j){ml, + ml+{m] - s) + (-m^ + s)t - ml,{s + t - u))) 
+fiRR{j)m{-'2rn% + 2m^i + i(m| + mj - s - u) 
+m^+(-m| - + s + u)))) + /iLL(?)(3m-+m^(/ii,R(j)m^o(m| 
-t) + 2fiRR{j)met) + fiRL{j)mem^o{-2ml, + 2m^t + t(m| + 
+m|+ (-m| - + s + li)) + fiLL{j){'^rn^rn% + m-t - mjt 
—2m^m'^t + m^t + m^si + rrv^st — rrv^t^ + rriftu — rn^tu + t^u 
+m|+ {—Tnj^ + mj — tu — m^{s + u) + mj^ {s + t + u)) — mj^ (2m^y 
-2m^t + t{s + t + u)))) + /iRij(i)(3m^+m^(/iRi,(j)m^o(m| - mj 
+'2fiLL{j)met) + fiLR{j)mem^o{-2ml, + 2m^i + i(m| + - s - 
+m^+ - + s + li)) + /iRR(j)(2m|m^ + m|^i - m^t 

— 2m^m^t + m^t + m^st + m^st — mlyt^ + m^tw — m\rtu + t^-u 
+m|+ {—TTij^ + mj — tu — ml{s + u) + (s + i + u)) — mj^ (2m^ 
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-2mlt + t{s + t + u)))))/{ml{mlo - i)(-m|o + t)) 
TiiXTii = ((4/2L/2Km^+m^ + /2l("^|+ - «) + - '"))("^| 

+ (m^ - uf - 2m|(m^ + u)))l{ml,{m\ - uf) 
TixTii = {2m^o{fiLR{i){f2Rmi{ml, - 2m^s + m|(3m^ - " 

+m^+ (m| + 3m^ - xi) - m^w + tu) + f2Lm^+{-rn^^ + m^s 
+m| (-m^ - 3m^ + s) + 2m^t + m^-u - s-u + m^(-3m^. + u,))) 
+fiRL{i){f2Lme{ml, - 2m^s + m|(3m^ - ^) - "^^^ + "^^+("^| + 3"^^ 
-li) - m^ii + iw) + /2iim-+ (-m^ + m^s + m| (-m| - 3m^ + s) 
+2m^t + m^M - + mji-Sml, + u)))) + fiLL{i){'2f2Rm^+me{mj^ 
+2m^ - m^s - 2m^t + m| (m^ + - s - 2u) + m^{3m^^, - u) 

o 2 I , 2\ , r /r- 2 2 I 4 22, 22j. 4j. 

—6m^u + su + u j + /2L(5m£m^yS + m^s - m^^s + m^m^t - m^t 

i2j.2i4/j. \ 22 4 2 2 ,2 

+m^r + m^^ (s — r — xij — m^m^u — rriy^u — m^su — rriy^su -\- s u 

—rnjtu + rn^tu — t^u + m^w^ + 2m'^u^ — + (— + 

+m^(— 8m^ + s + t — — sti + + 2ti^ + 2m^(s — t + u)) 

+m|+(2m|^ - (m^ - u){—s + t + u) + mj{—6m^ + 2u) 

-m|(2m| + 2m^ - s + i + Su)))) + /iRR(i)(2/2Lm^+m^(m|^ + 2m^ 

—m^^s — 2m^^t + (m^ + — s — 2-u) + mj{3ml^ — u) — 3m^M + su 

i2\iX/'r-22 I 4 22, 22j. 4j.i 2j.2 

+u ) + j2R{5mgm^s + m^s - m^s + m^m^t - m^t + m^t 
+m|^ {s — t — u) — m^rn^u — m^u — m^su — rn^su + s^u — m^tu 
+m^tu — t^u + m^M^ + 2m^M^ — + rnj^ {—s'^ + 
+mj{-8ml, + s + t- u) - su + tu + 2v^ + 2m^,(s - t + u)) 
+m?+(2m|^ - (m^ - u){-s + t + u) + m|(-6m^ + 2u) 
-m|(2m| + 2ml -s + t + 3ii)))))/(2m^(m|o - t){ml - u)) 
\T\' = Y.'^ixTi + TuxTu + 2j2TixTu (C9) 
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I. t- channel 234) exchange 

II. K-channel exchange 

Note the L-R switch for the couphngs below. 





- ^h-e- 




flLR{i) 




(^R 

_~o ^H+-xT-X° 


fiRhii) 


- ri 


_~o ^H+-xI-X° 


flRR{i) 


- 

ei-e- 


fyR 

_~o ^H+-xT-x'i 






-Xi ^Vi-h-H+ 






-xT ^vi-li-H+ 




= i^x^i 


— flRL{i){flRR 



+/iRL(i)(r«|+ u))) + fiLL{j){fiLR{i)m^+{-mj^ +m'j + t) 

+fiRL{i)m{'^rnl+ +m'~^+m'j - s- u)) - fiLR{i){m^o{2fiRL{j)m^+mi 
+fiLR{j){m^^+ + m| - u)) - fiRR{j)me{2m'i+ +mj^+rnj-s- u))) 
+fiLL{i)ifiLRij)m^+m^o{-mj^ + m'j + t) + me{-2fiRR{j)m^+t 
+fiRL{j)m^o{2mi+ + m| + - s - u)) + fiLL{j){mj^ - + nT-h 
+m?,+ (2m|^ — 2m| —t) + st + rriju — ?7i|^(s + t + u))) 
+fiRR{i){fiRL{j)m^+m^oi-m\ + + t) + mi{-2fiLL{j)m^+t 
+fiLR{j)m^o{'2ml+ +m1^+mg - s- u)) + fiRR{i){m\^ - mj + m^s 
+m^+(2m| - 2m^ -t) + st + mju - m|(s + i + u)))) 
/((m|o-i)(-m|o+i)) 
TJixT^i = {Af2Lf2Rm^+mi + flL{mi+ + m| - + fm{nvi+ + mj - u))/{ml^ - uf 
TixTii = {m^o{fiLR{i){2f2Rm^+rrn + /2L(m|+ + - u)) + fiRL{i){'2f2Lm^+m 
+/2ii(m?+ + - + fiRR{i}(f2Rm^+{mj^ - - t) 
+f2Lmi{-2mi+ - m| - m| + s + rt)) + fiLL{i){f2Lm^+{m'j^ - mj - t) 
+/2R"^K-2rn^+ - m| - + s + m))) 
/((m|o-i)(m?^-K)) 
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|r|2 = ^7ix7i + 7Jix7ii + 2^7ix7ii (CIO) 



Appendix D: Neutralino-Sneutrino Coannihilation 

Below is the list of the amplitudes squared for neutralino-sneutrino coannihilation. Note 
that, for identical-particle final states, one needs to divide them by two when performing the 
momentum integrations. Below, v refers to 

vv* — > WW 

I. s-channel H annihilation 

II. s-channel h annihilation 

III. M-channel exchange 

IV. point interaction 

V. s-channel Z annihilation 



/l 


— Ch-w-w C~_~_fj 




/2 


= Ch-w-w C~_~_h 




h 






h 


— ^v-v-W-W 




h 


— Cz-w-w C~_~_z 






= (12mt - 4m^s + s')/(4mt( 


s - mlf) 




= (12mt-4m^s + s')/(4mt( 


's-mlf) 




= {ml + {ml, - uf - 2m?(m^ 


+ u)fl{mi{u-mlf) 




= (12mt-4m^s + s2)/(4mt) 


1 




= (-32m|m^ - 24m^s + (i" - 


- v^f - 8mt(s^ - 2(i - uf - s{t + u)) 




+2ml{s'^ -2s^{t + u)-2{t 


- uf{t + u)- s{f - 6tu + u^)) 



+4m^(16mt + {t - uf + m^(-6s + 8{t + u))) 

+4m?(24m^ - {t - uf{t + u)+ 4mt(s - 2{t + u)) + 2m1{-Atu 

+s{t + um/i^<is-mir) 
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TJxTJi = (12mt-4m^s + s')/(4mt(s-ml,)(s-m^)) 

7^x7^11 = {Gm'l + miiGml, - s) - su"^ + 2ml,u{s + 2t + u) - ml^ibs + At + 8u) 

+2m?(10mt + SU- m^(s + 2t + Au))) / {2m'l{s - mjj){u - mjj) 
TjxTiy = -(12mt-4m^s + s')/(4mt(s-m^)) 

TixTy = -{{t - u) (16mt + mliAmI - 2s) + s{t + u)- 2m\{2s + t + u))) 

/(4mt(s-m^)(s-m|)) 
^ixTJn = (6m^ + m^(6m^ -s) -s'u^ + 2m^M(s + 2t + 'u) -mt(5s + 4i + 8M) 

+2m?(10mt + sii - m^(s + 2i + 4ii)))/(2mt(s - ml){u - mlj) 
TuxTjy = -{12mi-4mls + s^)/{4mi{s-ml)) 

TuxTy = -{{t-u){16m'^„ + ml{4m1-2s) + s{t + u)-2m1{2s + t + u))) 

/{im'lis - ml){s - ml)) 
TnixTiy = -(6m^ + m^(6m^ - - si*^ + 2m^ii(s + 2i + -u) - mt(5s + 4i + Sii) 

+2m2(10mt + - m2,(s + 2t + 4M)))/(2mt(M - m?J) 
TJiixTV = (tV-'u^ + 2m|(8m^-t + M) +4m^(s-t + 'u) + mt(2s2 + 3t^ 

+4s(i -u)+ Atu - 7u^) + 2ml^u{s^ - 2st - 2t^ + 2u^) - 2m?(8m^ 

+mt(4s + 5t- 5u) + u{t^ + tu- 2u^) + m^(s^ - 2st - 2t^ + 2su - Stu 

+2u^)) + mlif + Atu - 51? + 4m^(s - A{t + u)))) 

/{2ml,{s - ml){u - mlj) 
TjyxTy = {{t-u){16ml, + ml{Am'l-2s) + s{t + u) -2ml{2s + t + u))) 

/{4mi{s-ml)) 

= f^TjXT, + flTuXTu + f!TuiXTui + f!T,yXT,y + fiTyXTy + 2fj2TlXTu 

+2fif3TixTni + 2/i/47ix7iv + 2hfaiX% + 2f2h%iX%ii + 2f2hruxTiy 

+V2hr^ix-ry + 2fshTmxrw + Vzh%iix-ry + 2UU%Yxry (di) 

I. s-channel H annihilation 

II. s-channel h annihilation 

III. ii-channel v exchange 
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IV. i-channel v exchange 

V. point interaction 

/i = Ch-z-z C~_~_jj 

f2 — Ch-z-z C~-Z-h 

h = ic~_~_zy 

TixTi = {12ml- Amis + s'^)/{Aml{s-m]jf) 
TuxTii = {I2ml - Amis + s^) / {Aml{s - mlf) 
rmxTm = {mt + {ml-uf -2ml{ml + u))y{mi(u-mlf) 
TjyxTiy = {mt + {ml-tY -2ml{ml + t))y{ml{t-mlY) 
TyxTy = {12ml - Amis + s-^) / (Ami) 
TixTu = {I2ml - Amis + s^) / {Aml{s - m]j){s - ml)) 
TixTju = {6ml + mt{6ml - s) - su^ 

+2mlu{s + 2t + u) - ml{5s + At + 8u) + 2m?(10m^ + su - ml{s + 2t 
+Au)))/{2ml{s-ml){u-ml)) 
TixTiY = {6ml + mU6ml- s) - st'^ + 2mlt{s + 2u + t) -ml{5s + Au + 8t) 

+2m?(10mt + st- ml{s + 2u + At)))/{2ml{s - m]j){t - ml)) 
TixTy = -{12ml- Amis + s^)/{Aml{s-mjj)) 
TiixTiii = {6ml + ml{6ml - s) - su^ + 2mlu{s + 2t + u) - ml{5s + At + 8u) 

+2m?(10mt + SU- ml{s + 2t + Au))) /{2ml{s - ml){u - mj)) 
TijxTiy = {6ml + mt{6ml- s) - st^ + 2mlt{s + 2u + t) -ml{5s + Au + St) 

+2ml{10ml + st- ml{s + 2u + At)))/{2ml{s - ml){t - mj)) 
TiixTy = -{12ml- Amis + s'^)/{Aml{s- ml)) 
TuixTiy = {ml + ml + m?(6m| -t-u)+tu- m|(2s + t + u))'^ 

/ {ml{t - ml){u - ml)) 
TuixTy = -{6ml + ml{6ml - s) - su'^ + 2mlu{s + 2t + u) -ml{5s + At + 8u) 
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+2m|(10mt + su — m\{s -\-2t-\- Au)))/{2m2{u — m|)) 
TiyxTy = -{6ml + mi{Qml- s) - st"^ + 2mlt{s + 2u + t) -mi{5s + 4:U + 8t) 
+2m?(10mt + St - ml{s + 2u + 4t)))/(2mt(t - m?)) 
|r|2 = f^TixTj + f^TuxTu + fiTjuxTju + f!TiyxTjy + fiTyxTy + 2fj2rixTu 
+2hhT^xTni + 2hfJ{x%y + 2hh%xTy + 2hh%^x%n + 2hfaiiX%j 
+2/2/5^ixrv + 2fsUTuixTjy + 2hf^TuixTy + 2/4/57Jvxrv (D2) 

I. s-channel Z annihilation 

II. i-channel exchange 

III. w-channel exchange 

/2 = C~_~_f^C~_~_z [C~_~_fj C~_~_z] 

TixTi = {-16mlml + Amis + (t - m)^) /{Aml{s - mlf) 
T^xTn = (mt + (ml-ty -2ml{rnl + t))/(ml(t-miy) 
TuixTui = (m4 + (m|-'fi)^-2m?(m| + '?i))/(m|(?x-m?)^) 
T^xTJi = -{t{t -u)+ ml{-8ml - t + u) + rr4(2s - t + M))/(2rr4(s - ml){t - ml)) 
TixTui = -(u(u-t)+ml{-Sml-u + t) + ml{2s-u + t))/{2ml{s-ml){u-ml)) 
TiixTui = (-m^ + m^ + iii-m?(-6m^ + i + ii) -m^(2s + i + ii)) 
/(m|(?i - r74)(t-m?)) 
\r\' = f^T,xT, + f^TuxTu + fiTuixTui + 2fj2TixTu + 2fJ,TixTui 

+2/2/3^ix7i„ (D3) 

I. s-channel Z annihilation 

fl = C~_~_zCz-h-A [C~_~_zCz-h-a] 
%x% = (t-u)^/{s-ml)^ 

80 



ff 

I. s-channel Z annihilation 

fn = C~_~_z Cz-f-f 
%xT, = -{fl + fl){Amls-s^ + {t-uf)/{s-mlf 

For quarks |Tp is multiplied by 3 for color. 
Pi)* > tt 

I. s-channel Z annihilation 

II. s-channel h annihilation 

III. s-channel H annihilation 

flL = C~_~_z Cz-t-t 
flR — C~_~_z Cz-t-t 

fs — C~_~_fj Cn-t-t 

T.xTj = {^-2mlUiL-hR?s + Aml{2mlU\L-hR?-UlL + flR)s) 
+ {fiL + fiR){s'-{t-un/{s-mir 

T^xTn = 2(s - 4m2)/(s - m^)2 

TuixTui = 2(s-4m?)/(s-m^)' 

T.xTn = -{{hL + fiR)2m,{t-u))/{{s-ml){s-ml)) 

TixTui = -{{fiL + fiR)2mt{t-u))/{{s-ml){s-mj,)) 

T^xTui = 2(s-Amj)/((s-mjj)(s-ml)) 

|r|2 = 3(71x71 + f^TuxTn + fiTnxTn + 2h%xru + 2hTixTin 

+2/2/3^lX7i„) 
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vv* — > ee 



I. s-channel Z annihilation 

II. i-channel charginos exchange 



flR = C~_~_z Cz^e-e 



2 \2 



rix% = {-UtL + ftR){^<s-s'-r{t-uY))/{s-miy 

TuxTn = mf2{3){tu-ml)/{{t-m\^){t-ml+)) 

= {{l/2)h{3)hr^{Amls-s' + {t-uf))/{{s-ml){t-ml,)) 

|r|2 = 7ix7i + ^7iix7Ji + 2^7ix7ii (D7) 



z/z/ — > z/z/ 



I. s-channel Z annihilation 

II. i-channel neutralinos exchange 



/l — C^-v-Z ^Z-v-v 

/2(i) = {C~_~^_^f 

TJxTI = (-4m?s + s=^ - (t-M)2)/(s-m^)2 

TIixTIi = (-2m^ + 2tM)/((t-m|o)(t-m|o)) 

TixTJi = _(i/2)(-4m?s + s2-(i-ii)2)/((s-m|)(i-m|o)) 

|r|2 = f^%xTj + J2f2{i)f2{j)TuxTu + 2j2fif2{i)TixTu (D8) 



z/z/* — > 

I. s-channel H annihilation 

II. s-channel h annihilation 

III. ^-channel cl exchange 

fl — Ch-w+-h- C~_~_jj 

f2 = Ch^w+-H- C~_~_f^ 
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/3 — ^u-CL-W ^v-eL-H+ 

TixTi = (m^+ - 2m]j+ {2ml + + s-t-u) + (-2m? + - s + i + uf) 
/{Aml{s-mlf) 

TnxTu = {m%+ - 2m%+ (2m? + 7m^ + s-t-u) + {-2m\ + m^ - s + i + uf) 
/{Aml{s-mlf) 
TnixTni = [m^ + {ml - tf - 2ml{rnl + t)) / {ml{t - ml^f) 
TixTu = {m%+ - 2m%+ (2m? + 7m^ + s-t-u) + (-2m? + m^ - s + i + -uj^) 

/(4m^(s-m^)(s-m^)) 
TixTui = (2m^ + mt - m^s - 2m1t -st + t^ + m^+ {-mj - 3m^ + t) 
+ml{m1 + s - 3t - u) + m^u + tu)/{2ml{s - m\){t - mjj) 
%iX%ii = (2m^ + - m^s - 2m^i - si + + m^+ (-m? - 3m^ + t) 
+ml{ml + s — 3t — u) + m^u + tu)/{2m1{s — ml){t — m\^)) 
|T|2 = f^TixT, + f^Tnxru + fiTuixTui + 2fJ,T,xTu + 2f,hTixTui 

+2f2h%ix%ii (D9) 

z>z>* — > H+H- 

I. s-channel H annihilation 

II. s-channel h annihilation 

III. i-channel exchange 

IV. point interaction 

V. s-channel Z annihilation 



/l 


— Ch-h+-h- 


/2 


= Ch-H+-H- Cu-u-h 


/3 




h 


— ^v-v-H+-H- 


h 


= C~_~_fjCz-H+-H- 


%x% 


- y{s-mir 


TuxTu 
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TmxTm = l/(i - m?J' 

TJvxTJv = 1 

TVxTV = (M-t)7(s-m2)2 

TixTii = l/((s-m^)(s-m^)) 

TJxTin = l/((s-m2^)(t-m?J) 

T^xTJv = -l/(s-m2,) 

TixTv = -{u-t)/{{s-m]j){s-ml)) 

ru-xTni = l/((s-m^)(i-m?J) 

TJixTJv = -l/(s-m2) 

TJixTv = -(«-t)/((s-m^)(5-m2)) 

TiiixTJv = -l/(i-m?J 

T^iixTv = -{u-t)/{{t-ml^){s-ml)) 

TJvxTv = («-t)/(s-m2) 

|r|^ - /i'7Ix7i + /27iix7ii + /|7i„x7i„ + /27Ivx7iv + /5Tvxrv + 2/i/27ix7ii 
+2fJsTixTni + 2/i/47Jx7Jv + ^fJ^ixTy + 2f2h%iX%n + ^hfaiixr^ 
+2/2/5^ixrv + 2hfam~x%j + 2fsfauixry + 2Ufajyy~% (DIO) 

^ if if [/i/i] [/iif] 

I. s-channel if annihilation 

II. s-channel h annihilation 

III. t-channel V exchange 

IV. point interaction 



/l 


— Ch-h-h C~_~ 


-H 




-h-h C~_~_h] [Cn-H-h C~_~_ 


/2 


= Cn-H-h C~_~_ 


-h 


[Ch-h- 


-h C~_~_f] [Cn-h-h C~_~_,^ 


/3 


= {Cz-Z-h)'^ 




-u-h)'^] 


[Cp-v-h Cu-u-h] 


U 


— C~_~_H-H 




ii—h—h\ 


[Cu-u-H-hl 




= l/(«-m|r 








TuxTu 
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i/(i - mir 




1 






!/((«- 


m\){s — m\ 


TixTui = 


!/((«- 


m]j){t-rn\ 






-m]j) 


TiixTiii = 


l/((^- 






-l/{s- 


-ml) 


^iixTJv = 


-1/it - 


-ml) 


|T|2 = 







+2fJsTixTni + 2f,UTixTiy + 2hhTuxTni + 2hhTiixTv, 
+2/3/4^iix7iv (Dll) 



uu* — > AA 



I. s-channel H annihilation 

II. s-channel h annihilation 

III. point interaction 

/i = Ch-a-a 

/2 = C~_~_f^ Ch-A-A 

%x% = l/{s-m]jf 

TuxTii = l/{s-mlf 

TnixTui = 1 

%xTn = l/{{s-mj,){s-ml)) 

TixTui = -l/(s-ml) 

TuxTni = -1/is-ml) 

|r|2 = flT,x% + flTuxTu + fiTuixTui + 2fj2TixTu + 2Ma,xTui 

+2f2fsTuxTjn (D12) 
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AZ 



I. s-channel H annihilation 

II. s-channel h annihilation 



/i 



C~_~_^ Ch-Z-A 

{m\ - 2m\{2ml + 7ml + s - t - u) + (-2m? +ml-s + t + uf) 
/ {Aml{s - m]jY) 

(m\ - 2m\{2ml + 7ml + s - t - u) + (-2m? + ml-s + t + uf) 
/{Aml{s-mlf) 

{m\ - 2ml{2ml + 7m| + s-t-u) + (-2m? +ml- s + t + uf) 

/ {Aml{s - m]j){s - ml)) 

fl%x% + flTuxTu + 2/i/27ix7ii (D13) 



vv 



I. i-channel neutralino exchange 

II. li-channel neutralino exchange 



m 



= (2sm^om^o)/((t - m|o)(t - m|o)) 

= (2sm^omj^o)/((M — m|o)(M — m|o)) 

= (2sm^om^o)/((i - m|o)('u - m|g)) 

= E/W/(i)(^x7J + 7Jix7Ji + 27ix7ii) 



vZ 



I. s-channel v annihilation 

II. i-channel V exchange 

III. w-channel neutralino exchange 



(D14) 



/i 



C'~_~o_^ Cz-v-v 



86 



TjxTi = {l/2){l/ml){mp-mls + s{-ml-st + ml{s + t-u)) 

-m\{2m\ - 2mls + s{t + it)) + ml{2m\ - 2m|s + s(s + i + u))) /{sf 
TiixTii = {l/2){ml-t){mt + {ml-tf -2ml{ml + t))/{ml{t-mlf) 
TuixTju = (l/2)(-l)(m^o/3i(i)(3(-m^)m2/3s(j)(-m| - ml + s + t) 
+'m^of3L{j){ml + s(-m| + u) - ml{s - t + u))) 
+f3R{i){-H-ml)m^omlf3L{j) - m|/3ij(j) 

+3{-m^)m^omlfsLij){-ml + s + t) + m'^fsji{j){-2ml + s + t + 2u) 

+/3i?(j)(2m^ + 2mlsu + tu^ - 2m^(s + i + u)) 

+mlf3R(j)(mi + ml(s + t + u) - u(s + 2t + u)))) 

/{ml{u - m\o){u - m|o)) 

TixTii = {l/2){-l/{2ml)){mt{s + t-u) + {ml-t){-s'^ -t'^ + ml{s + t-u) + u'^) 

+m?(-s^ -st- 2f + 2ml(s - t - u) + tu + u^) + m|(m|(-s + t-5u) 

+ml{8ml -s + t-u)+t{s-t + u))) /{s{t - ml)) 

TixTni = {l/2){l/{2ml)){-2{-ml)m^ofsL{j)s-mysR{j){s + t-u) 

+2{-my^)m^of3L{j){2ml - 2m|s + s(-m~ + t + u)) 

+mlf3R(j)(s^ + 2st + + su + tu- 2u^ + m?(8m| - 3s - i + tt) 

-ml(5s + t + 3u)) + /3K(j)(-4m^s + mjiAml - 4ml{s + t) + {s + t- u)u) 

-u{s^ + - v^) + m|(4s2 + it- u)u + s{At + u))))/{s{u - m|o)) 

Xj 

%ix%u = (l/2)(l/m|)(mj/3ij(j)(m? + 3m| - i) + (-mj^)m^o/3L(j)(m^ + m?(3m| - s) 
+si — m|(s + i + 2ii)) + ?77.|/3ij(j)(3m^ + m|(5m| — s — i — w) 

+ t + m) - m|(4s + 5t + 2ix)) + fm{j){-ml{s + t-u)-t^u 
+ml{s^ + 2st + f + tu- u^) + m?(tM + m|(-s - 3t + «)))) 
/((i-m?)(it-m|o)) 
|r|2 = f^TjxTj + f^%jxTu + J2'ruixTui + 2fj2TjxTu 
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+2 J2 (/i^xTJii + f2%ix%n) (D15) 

j 

Xi^ — > eW^ 

I. s-channel v annihilation 

II. i-channel e exchange 

III. li-channel chargino exchange 



/l 


= a- 


-e- 






= c~- 


-e- 


-W'-'ei-x°-e 


f2R{i) 


= c~_ 


-e- 




fsLii) 








/3r(0 






"-e^X?-X+-W^ 



TixTi = (l/2)(m|s - m^s + s(-m^ - + m^(s + t - n)) - m|(2m^ - 2m^s 
+ u)) + m?(2mt - 2m^s + s{s + t + u)))/{ml,s^) 
TuxTu = (l/2)((/2i^(i)/2«(i) + /2LW/2L(j))(m|-i)(m^ + (m2^-if 
-2ml{ml + t)))/{ml{t - m?J(i - m|)) 
^iixTIii = (-l/2)((/3L(2)m-+(3(-m^)m'^/3i?(j)(-^ - ml + s + t) 
+rn^+ hL{3){rn^w + ^{-m\ + u) - m1{s - t + u))) 
+f3R(i)(-^-m^)^m~+m1f3L(j) - m\fm{3) 

+S{-m^)m^+ml,f3Lij)i-m1 + s + t) + m|/3/j(j)(-2rr4 + 5 + ^ + 2^) 
+/3R(j)(2mt + 2mlsu + tu^ - 2mt(s + t + uj) 
+m|/3i?(j)(mt + mlis + t + u) - u{s + 2t + u))))) 
/{ml{u - ml+){u - m|+)) 
TixTu = {-l/2)U2R{j){mUs + t-u) + {ml-t){-s^ -t^ + mlis + t-u)+u^) 
+ml{-s'^ -st- 2t^ + 2m^(s - t - u) + tu + u^) + m\{ml{-s + t-bu) 
+ml{8ml, - s + t-u) + t{s-t + u)))) / {2m1s{t - ml.)) 
TixTJii = {i/2){{-2{-m^fm^+hL{j)s-m\hR{j){s + t-u) 
+2{-m^)m-+f3LU){2mi - 2m^s + s(-m? + t + «)) 
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+mlf3R{j){s^ + 2st + f + su + tu- + m\{%m\, - 3s - i + 
-m^(5s + t + 3m)) + fprj{-4m^s + ml{Aml, - Aw?^{s + t) 
+ (s + t - u)u) - u{s^ + - V?) + m^(4s^ + - «)« + s{At + «))))) 
/{2ml,s{u - m|+)) 

^ixTiii = (l/2)(/2B(i)(mJ/3R(i)(m| + 3m^ - i) + (-m^)m-+/3i:(i)(mt + mi(3m^ 
-s) + - m^,,(s + t + 2m)) + m|/3i?(j)(3mf,, + m|(5m^, - s - t - m) 

+ i + m) - m^(4s + 5i + 2u)) + f3R{j){-mt^{s + t-u)-fu 
+m^(s^ + 2si + + - m^) + ml{tu + m^(-s - 3i + m))))) 
/(m^(t-m|)(M-m^+)) 
|r |2 = /^TixTI + ^ TiixTii + ^ TiiixTin + 2/i ^xTJi + 2/i ^x7I„ 

+2^7iix7J„ (D16) 



— > hv [Hv] 

I. f-channel neutralino exchange 

II. M-channel u exchange 

= C'~_~o_^ Cjo_~o_ft [C~_~o_^C^o_~o_h\ 

/2 = Cx;-j^o_^ C~_~_ft [C~_~o_^ C~__~_^] 
TixTi = {l/2){fiR{i)m^o{fiR{j)m^o{ml - u) + hL{3){-m^){-m\ - ml^H] + s + u)) 
+fiL{i){fiR{j){-m^)m^o{-ml - m^H] +s + u) 
+/il(j)(^ - mh[H\ - St- m|(s + u)+ ml^H]{s + t + u)))) 
/((t-m|o)(t-m|o)) 
TiixTu = {l/2){ml-u)/{u-mlf 

TixTii = (l/2)(/iR(i)m^o(m| - u) + fiL{i){-m^){-ml - m\jj] + s ^ u)) 

/{{t-mlo){u-ml)) 

i,3 i 
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I. t-channel neutralino exchange 



fiii) - C~_~o^^C~o_~o_A 

TixTi = (1/2) (/^i(i)m^o(/^i(j)mx°(^-'") + /iO')(-^x)(-^-^i + « + 

+/L(i)(/R(i)(-m^)m^o(-m| -ml + s + u) + fL{j){m\ -m\-st- m\{s + u) 
+ml{s + t + - m|o)(t - m|o)) 

|r|2 = (D18) 



2J 



I. t-channel chargino exchange 

II. xi-channel cl exchange 



haR{i) = C'|_j-_^+ 

fibR{i) = C~_~+ 

/2 — C~_~o_^ C~_~_jj+ 



%x% = (l/2)(/i6i?(i)/i6R(j)K/i„i(i)/i„i(j)-m|,+ /i„i(i)/i„i(j) 

\3/ 



-(-m^^) {m^+flaL{i)flaR{j) + m^+flaL{j)flaR{i)) " flaL{l)flaL{j)st 

-{-'m'x){m^+fiaL{i)fiaR{3) + m^+fiaLij)fiaR{i)){mH+ -s-u) 
-m~+m~+fiaR{j)fiaR{i)u + m]j+fiaL{i)haL{j){s + t + u) 

+ml{m~+m~+fiaR{j)flaR(i) - flaL{i)flaL(j)(s + u)))) / ((t - - 1711+)) 

TuxTu = (l/2)(m|-M)/(M-m2j2 

TixTii = {l/2){fibR{i){-{-m^)^fiaL{i)+mlm^+fiaRii)-m^+fiaR{i)u 
+i-m^)fiaL{i){-m]j+ + s + u)))/{{t - m|+)(xi - m?J) 
|r|2 = ^TjX% + f^TuxTu + 2f2^TixTu (D19) 
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I. t- channel neutralino exchange 

f{i) = {C~-^o_^) 



%x% = f{i)f{j){-2mt + 2tu)/{{t-mlo){t-mlo)) 
|r|2 = ^T^xTJ (D20) 



I. t-channel charginos exchange 



%x% = f{i)f{j){tu-ml)/{{t-ml,){t-m\,)) 

i,3 



I. f-channel neutrahno exchange 



i2 



TJxTJ = /(i)/(j)(2sm^om^o)/((f-m|o)(i-m|o)) 
|r|2 = TJxTi (D22) 
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